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CONVEXITY  OF  DOMAIN  FUNCTIONALS 


by 

P.  R.  Garabedian  and  M.  Schiffer 

CHAPTER  I 
INTRODUCTION 


A problem  of  considerable  significance  is  the  study  of  the  dependence 
of  the  Green's  function,  Neumann's  function,  and  eigenfunctions  of  a linear 
elliptic  partial  differential  equation  on  their  domains  of  definition.  The 
importance  of  this  question  lies  in  the  difficulty  which  is  generally 
encountered  in  the  explicit  calculation  of  such  functions.  Since  elementary 
formulas  are  the  exception  rather  than  the  rule  In  the  theory  of  partial 
differential  equations  in  an  arbitrary  domain,  one  turns  to  the  investigation 
of  the  properties  of  the  basic  solutions,  and  variational  formulas  exhibiting 
the  domain  dependence  of  these  solutions  furnish  one  of  the  principal  tools 
of  the  investigation. 

One  of  the  most  common  applications  of  the  theory  of  variation  of  the 
Green's  function,  Neumann's  function  and  eigenfunctions  arises  in  the  study 
of  extremal  problems  for  the  capacity,  virtual  mass  and  eigenvalues  of  a domain. 
We  are  led  to  such  extremal  problems  on  the  one  hand  in  an  attempt  to  estimate 
these  domain  functionals  in  terms  of  the  geometry  of  the  domain,  and  on  tne 
other  hand  by  the  equivalence  of  the  solution  of  particular  variational 
problems  with  the  existence  and  uniqueness  of  solutions  of  physical  problems, 
such  as  the  construction  of  free  boundary  flows.  Extensive  investigations  in 


these  two  directions  have  been  carried  through  for  Laplace's  equation  in  the 
plane  and  also  for  other  typical  elliptic  partial  differential  equations  in 
two  independent  variables  [2,  7,  8,  9,  18]. 

In  this  paper,  we  develop  a rigorous  theory  of  variation  of  domain 
functions  in  space  of  three  dimensions  as  well  as  in  the  plane.  We  not  only 
present  an  adequate  mathematical  discussion  of  the  classical  Hadamard 
variational  formulas  in  space,  but  also  generalize  the  so-called  interior 
variational  method  to  three  dimensions.  We  derive  expressions  for  the  second 
variations  of  the  capacity,  virtual  mass,  and  other  physical  quantities 
already  mentioned,  and  we  deduce  from  them  various  interesting  convexity 
theorems  for  these  domain  functionals . 

Our  investigation  of  the  second  variation  was  motivated  by  the  suggestion 
of  Max  Shiffman  that  in  cases  where  one  can  guess  a domain  for  which  a certain 
combination  of  domain  functionals  is  stationary,  one  might  be  able  to  apply 
a minimax  theory  in  order  to  prove  that  the  domain  in  question  actually 
minimizes  that  combination,  provided  one  could  show  that  whenever  the  first 
variation  of  the  combination  vanishes  the  second  variation  is  positive-definite. 
This  suggestion  i3  based  on  the  fact  that  on  certain  surfaces,  corresponding  to 
two  points  whose  heights  are  relative  minima  there  exists  a saddle- point . While 
we  have  not  had  any  direct  success  with  this  line  of  reasoning,  we  have  been 
able  to  deduce  a number  of  uniqueness  theorems  from  convexity  properties  of  the 
domain  functionals  which  are  based  on  the  second  variation.  One  should  mention 
in  this  connection  the  work  of  Friedrichs  T 6 ],  who  proved  by  such  a method 
the  uniqueness  of  certain  free  boundary  flows. 

In  Chapter  II  we  define  interior  variations  of  a 3 -dimensional  domain  D 
by  means  of  differentiable  mappings  cf  D depending  on  a small  parameter  £ „ 
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The  first  order  shifts  in  terras  of  £ of  the  Green's  function,  Neumann's  function 
and  eigenvalues  which  result  from  this  variation  of  D are  calculated  rigorously 
by  referring  all  varied  Quantities  back  to  the  original  domain  D through  the 
infinitesimal  mappings.  The  study  of  the  domain  dependence  of  these  basic 
functions  for  a given  linear  elliptic  partial  differential  equation  is  thus 
transformed  into  an  investigation  in  the  fixed  domain  D of  the  dependence  of 
the  solution  on  the  coefficients  of  the  equation.  Such  an  investigation  is 
readily  carried  through  by  Hilbert's  classical  method  based  on  integral 
equations  [ill.  The  variational  formulas  which  result,  for  example,  in  the 
case  of  the  capacity,  are  given  in  terms  of  domain  integrals  involving  the 
Maxwell  tensor  of  the  relevant  electrostatic  field,  and  their  validity  does 
not  depend  upon  strong  smoothness  assumptions  on  the  boundary  of  the  domain. 
However,  when  the  boundary  of  the  domain  is  sufficiently  smooth,  our  interior 
variational  formulas  yield  by  application  of  the  divergence  theorem  the  classical 
Hadamard  variational  formulas,  for  which  we  thus  obtain  a strict  derivation. 

Our  3 -dimensional  variational  theory  collapses  easily  by  specialization  to 
the  better  known  theory  of  variation  in  the  plane. 

Once  in  possession  of  a rigorous  proof  by  the  Hilbert  method  that  we  can 
expand  the  varied  domain  functions  in  powers  of  £ , we  are  justified  in 
employing  the  perturbation  method  to  calculate  the  second  variations  of 
these  functions.  We  do  this  in  Chapter  III  and  obtain  interesting  second 
variation  expressions  for  the  capacity,  virtual  mass  and  eigenvalues  corresponding 
to  various  particular  ways  in  which  vie  can  shift  D.  A number  of  convexity 
theorems  for  these  domain  functionals  are  the  outgrowth  of  this  investigation. 

For  example,  if  we  shift  the  surface  of  D along  level  surfaces  of  a harmonic 
function  U,  the  capacity  of  B with  respect  to  a fixed  interior  point  turns  out 
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to  be  a convex  function  of  U.  Similarly,  the  capacity  of  a convex  domain  D 
is  a convex  function  of  a parameter  which  gives  a variation  of  D defined  in 
terms  of  the  support  function  of  the  surface  of  D.  These  convexity  theorems 
are  applied  to  establish  extremal  properties  of  domains  for  which  it  Is  known 
that  the  capacity  is  stationary  under  certain  constrained  variations,  and  thus 
a number  of  uniqueness  theorems  are  deduced c Such  a theorem  can  be  obtained  al30 
for  the  problem  due  to  Evans  [ 5 3 of  finding  a surface  of  least  capacity 
enclos ing  a given  curve . 

In  Chapter  IV  we  specialize  our  variational  theory  to  the  case  of  two 
independent  variables  in  order  to  apply  it  to  show  the  existence  of  vortex 
sheets  in  axially  symmetric,  irrotational  flow  of  an  incompressible  fluid. 

An  indication  of  these  results  will  be  given  in  the  next  section,  where  we 
sketch  heuristically  an  extremal  characterization  of  vortex  sheets  in  3 -dimen- 
sional space  without  symmetry  of  any  kind.  In  the  axially  symmetric  case,  the 
convexity  of  the  virtual  mass  in  dependence  on  the  domain  can  also  be  used  to 
discuss  the  extremal  characterization  of  a vortex  sheet. 

Chapter  V Is  devoted  to  the  study  of  the  eigenfunctions  and  eigenvalues 

of  the  vibrating  membrane . Using  the  second  variation,  we  show  that  under 

certain  conformal  mappings  of  a domain  depending  on  a suitable  parameter,  the 

inverse  square  of  the  principal  frequency  of  the  domain  becomes  a convex 

function,  of  the  parameter.  This  theorem  is  applied  to  show  that  for  fixed 

principal  frequency  of  a domain,  the  capacity  with  respect  to  an  interior 

point  is  a maximum  when  the  domain  is  a circle  about  the  point,  a theorem  due 

n 

to  Polya  and  Szego  [l6]„  . 
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2.  Vortex  sheets. 

In  order  to  motivate  the  detailed  study  of  variational  methods  in  3- 
dimensional  space  which  follows  in  later  chapters,  we  treat  heuristically  in 
this  section  a 3 -dimensional  extremal  problem  for  the  virtual  mass  of  a 
steady  irrotational  flow  of  an  incompressible  fluid  whose  solution  yielas 
a flow  past  a vortex  sheet.  We  have  only  succeeded  in  proving  the  existence 
and  uniqueness  of  vortex  sheets  from  this  extremal  characterization  in  the 
cases  of  plane  and  axially  symmetric  flow,  but  a formal  discussion  of  the 
general  3 -dimensional  case,  which  is  the  one  of  interest  in  aerodynamic 
theory,  should  indicate  the  direction  in  which  further  development  of  the 
material  in  this  paper  should  be  pushed. 

Let  C be  a simple  closed  curve  in  space  and  let  £ be  a closed  surface 
which  intersects  every  curve  looped  around  G and  bounds  an  infinite  domain  D, 
Let  Cp  = <p(x,y,z)  be  a harmonic  function  in  D with  the  expansion 

(1.2.1)  Cp  = x + + . , . , 9 x2  + y2  + z2  , 

r 

near  infinity  and  with  an  inner  normal  derivative  which  vanishes  on  21, 

(1.2.2)  = 0 . 

dV 

The  function  (p  represents  the  velocity  potential  of  a steady  3 -dimensional 
flow  past  21  of  an  incompressible  fluid  in  D.  The  coefficient  &\  in  (2.l) 
is  related  to  the  kinetic  energy  of  the  flow  and  wc  shall  call  it  the  virtual 
mas3  of  with  respect  to  the  x direction. 

We  show  formally  that  if,  for  a fixed  curve  C,  the  surface  2T  is  so  chosen 

that 

(1.2.3)  9 minimum  , 


then  2H  reduces  to  a vortex  sheet  enclosed  by  the  curve 
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Let  be  another  surface  spanned  through  C whose  normal  displacement 

‘bis  from  ZT,  measured  from  any  point  of  ZT,  is  infinitesimally  small.  We 

* * 
denote  by  D the  infinite  domain  bounded  by  2L*  , we  denote  by  (2>  the  velocity 

* * 

potential  of  the  flow  of  the  form  (2.l)  past  ]EL  , and  we  denote  by  pf  the 

* 

virtual  mass  of  this  flow.  We  attempt  to  estimate  the  difference  Of  - Pf 
in  terms  of  the  normal  shift  S IS. 

If  R denotes  the  surface  of  a large  sphere  enclosing  2 and  2_*  , we 
find  by  Green's  theorem  that 


(1.2.4) 


-q{  - 


4TT 


Ilf 


a cp 


CD 

' 3 1/ 


- - <? 
t 


d CT 


since  the  integral  is  independent  of  the  radius  of  the  sphere  R.  Here  dCT- 
denotes  the  surface  element.  A further  application  of  Green's  theorem  yields 

(1.2.5) 


(A* - <A 


4TT 


Jjtvcpfdt-^^tvcp'-vq)]2  ar  , 


D-D 


where  d't  is  the  volume  element.  From  (2.5),  it  is  clear  that  the  first  order 
term  SrS  of  the  difference  CX  ~ CX,  considered  as  a functional  depending  on  'bV . 
is 

(1.2.6)  ^ ] \ (VCp^S^dcr- 

zr 

Let  us  non  suppose  that  2 is  an  extremal  surface  for  the  minimum  problem 
(2.3).  We  first  conclude  that  2 must  reduce  to  a single  sheet  containing  no 
interior  points,  since  according  to  (2.5)  the  virtual  mass  Q\  decreases 

2 

monotonically  as  the  flow  region  D expands.  Secondly,  we  notice  that  (V  Cp  ) 
must  have  identical  values  on  both  sides  of  the  extremal  sheet  2,  since  the 
first  variation  (2.6)  must  vanish  for  every  normal  shift  of  the  extremal  surface, 

r-/  r* 

and  since  a normal  shift  of  the  sheet  2 corresponds  to  values  of  bib  which 
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differ  in  sign  only  on  the  opposite  sides  of  TT.  Thus  in  the  extremal  case, 
(VCj>)2  is  continuous  across  'll,  which  implies  by  Bernoulli's  law, 

(1.2.7)  \ (VCj>)2+p  - const.  , 


that  the  pressure  distribution  p is  continuous  through  ZT*  But  the  continuity 
of  the  pressure  p is  precisely  the  physical  condition  which  characterizes  ZH 
as  a vortex  sheet.  In  space,  the  velocities  on  the  two  sides  of  a vortex 
sheet  are  permitted  to  have,  within  the  tangent  plane  of  the  sheet,  quite 
unrelated  directions,  but  their  magnitudes  must  be  equal.  In  3 -dimensional 
aerodynamic  theory,  such  sheets  are  introduced  to  account  for  the  lift  and 
drag  produced  on  an  airfoil.  We  have  succeeded  here  in  connecting  this 
concept  with  the  minimum  energy  principle  (2.3). 
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CHAPTER  II 

THE  METHOD  OF  INTERIOR  VARIATION 


x.  Sangealiiiaa- 

We  is  ant  to  develop  In  this  chapter  the  theory  of  the  variation  for 
Green's  functions  -with  respect  to  linear  partial  differential  equations  of 
elliptic  type  under  a variation  of  their  domain  of  definition.  For  the  sake 
of  simplicity  we  shall  treat  the  case  of  three  Independent  variables  x^ 

(i*  1,2,3)  which  vary  in  a domain  of  three-dimensional  space.  Let 
p(x^)  be  a continuously  differentiable  function  in  this  basic  domain  Dq 
and  consider  the  partial  differential  equation  of  elliptic  type 

(2„l.l)  L[u]  = V2u-  pu  » 0 , V*”  ■ 

Let  D be  a subdomain  of  D and  let  us  assume,  at  first,  that  D has  a 
smooth  boundary  C„  We  assume  further  that  the  only  solution  of  ( 1 » 1 ) 
in  D which  vanishes  on  C is  the  trivial  solution  u=0.  In  this  case,  there 
exists  a Green's  function  of  the  domain  D with  respect  to  (l.l).  This  Green's 
function  u(P,Q)  (with  F <3  ^ is  characterized  by  the  following 

three  requirements : 

a)  For  fixed  D,  G(P,Q)  is  a solution  of  (l.l)  as  a function  of  P 
and  twice  continuously  differentiable  in  D.  except  at  the  point  Q„ 

b)  The  function  G(P,Q)  - is  continuous  in  D,  has  uniformly 

bounded  first  derivatives  with  respect  to  P,  except  at  P=Q,  and  its  second 

derivatives  grow  at  most  like  r(P,Q)  ^ if  P approaches  Q,  where 

3 

(2.1.2)  r(P»Q)2  « YZ  (x.-  ?j2 

i-1  x U 

c)  G(P,Q)  vanishes  for  P€-C,  QeD. 


t 


l ZK 


_ o „ 


It  is  well  known  that  these  requirement;,  determine  the  Green's  function 
in  a unique  way  and  that  the  Green's  function  satisfies  the  symmetry  condition 

(2.1.3)  G(P,Q)  - G(Q,P) 

Thus,  in  particular,  Green's  function  is  a solution  of  (l.l)  in  dependence  on 
the  parameter  point  Q also. 

If  u(P)  is  an  arbitrary  solution  of  (l.l)  in  D and  is  continuous  in  D+C, 
it  can  be  expressed  in  terms  of  its  boundary  values  on  C by  means  of  the 
Green's  function  in  the  form 

(2.1.4)  u(P)  "if  "<Q)40' Q ’ 

n ° 


where  7/^  denotes  the  interior  normal  at  the  point  Q with  respect  to  the 
surface  C and  where  dC7“Q  is  the  surface  element  at  Q. 


We  can  also  solve  the  inhomogeneous  differential  equation 
(2-1.5)  ^u  - pu  ■ f (x^)  , u ■ 0 on  C , 

if  f(x^)  is  Holder  continuous  in  D+C.  In  fact,  the  solution  u(P)  can  be 
represented  in  the  form 


(2  ..1.6) 


rrr 

u(P)  = -j  ( 

l/i Jm 


G(P,Q)f(Q)d'C, 


where  dT^  is  the  volume  element  at  the  point  Q. 

This  result  leads  to  an  interesting  interpretation  of  the  quadratic  form 


(2ol,r7) 


rcf,fi 


G(P,Q)f(P)f(Q)d'Cpd'tQ 


D D 


In  fact,  we  obtain  from  (2.1.5/  and  (2.1.6)  by  means  of  Green's  identity 

(2.1.8)  rtf. 


r,f]  "JjJ^tVu)2  + pu2]dT: 


The  right-hand  3ide  is  just  the  Dirichlet  integral  of  u with  respect  to  the 
differential  equation  (l.l).  In  the  special  case  that  p(x,)  is  non-negative 
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in  D,  the  Green's  function  G(P,Q)  is  the  kernel  of  a positive- definite  form. 
This  fact  makes  the  theory  of  ( 1 - 1 ) particularly  simple  for  the  case  of  a 
non-negative  coefficient  p(P). 


2. 


We  consider  a three  times  continuously  differentiable  vector  field  in  Dq 

(2.2.1)  Si  = Si(xj)  i,j-  1,2,3, 

and  the  transformation 

(2.2.2)  xi  = xi(x^;C)  “X^+ES^x^)  > 


which  depends  on  the  real  parameter  £ . Given  a proper  subdomain  DC  D^,  we 

may  choose  a bound  b(D,D  ) > 0 such  that  for  I b(D,D  ) the  domain  D is 

0 o 

* 

mapped  topologically  onto  another  proper  subdomain  D c Dq  with  smooth  boundary 
surface  C . At  the  same  time,  we  can  choose  b(D,DQ)  so  small  that  all  domains 
D obtained  still  possess  a Green's  function  G (P,Q)  with  respect  to  the  differ- 
ential  equation  (l.l).  Our  principal  aim  in  this  chapter  is  to  express  G (P,Q) 
in  terms  of  the  original  Green's  function  G(P,Q)  and  the  transformation  vector 
field  (2.1). 

For  this  purpose,  we  refer  the  function  G (P,Q;  back  to  the  original 
domain  D.  Let  P (P)~x  (x  ),  Q (Q)=-  Cr".)  and  consider  the  function 

(2.2.3)  g(P,Qj  6)  » G*[P* (P),Q* (Q)] 


which  is  well-defined  in  D,  twice  continuously  differentiable  except  for  P«Q, 
and  which  vanishes  if  either  argument  point  lies  in  the  boundary  surface  C of  D, 

it 

Let  u (P)  be  an  arbitrary  solution  of  (l.l)  in  D and  let 


(2.2.4) 


* , * 


u(P; £ ) - u (P  (P)) 


be  the  corresponding  function  in  D.  We  have 

■*» 

9xm 


* 4 — A O v 

0x 


3 


i - 1 ^Xi  3x* 


(2.2.5) 


9 


J- 1,2,3 
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* , * 


The  function  u (P  ) yields  a stationary  value  to  the  Dirichlet  integral 


(2.2.6)  UVu  ) * pu  ]a-!T*  * J]  XI  ( A)2* 

d*  d U • 1 ivt  ®*i  3x* 


pn  r Q dz 


* * 
for  all  functions  in  D with  same  boundary  values  on  C . Here 

e^Xi^  " d(x^,3t^,xy 


(2.2.7) 


is  the  Jacobian  of  the  transformation  (2.2).  Consequently,  the  function 
u(P; £ ) must  satisfy  the  Euler- La grange  equation  for  the  right-hand  integral 


in  (2.6). 

We  introduce  the  notation 

(2.2.8) 

^xi 

Aik  " 0(xi}  — * 

3-1 

(2.2.8’) 

p(xi)  - e(xi)p(xjL) 

"k 


^ki 


and  we  obtain  the  following  transformed  differential  equation  for  u(Pj £ ) s 

3 


(2.2.9) 


'w\5^u*  ■*?-/“ 


CA.,  /Ou  » 0 


This  equation  is  satisfied,  in  particular,  by  g(F,Qj  e ) in  dependence  on  P 
for  P6D,  F/ Q. 

The  differential  equation  (2.2.9)  was  obtained  as  the  Euler- Lagrange 
equation  with  respect  to  the  Dirichlet  integral 

(2.2.10)  Q (u)  -|I(  XT  * fi 


l.k-1  lk  axi  8xk 
This  integral  reduces  obviously  to 

(2.2.10’)  Qq[u]  »Jjj[(Vu)2^- pu2ldt: 


in  the  case  £ * 0;  Q^Cu]  will  play  an  important  role  in  the  subsequent 


considerations . 
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The  preceding  transformation  of  (l.l)  into  (2.9)  permits  us  to  consider 
the  Green’s  functions  of  varying  domains  with  respect  to  the  fixed  differential 
equation  (l.l)  also  as  the  Green’s  functions  of  the  fixed  domain  D with  respect 
to  the  varying  differential  equation  (2.9).  In  this  way,  the  dependence  of 
the  Green’s  function  on  the  parameter  Z can  be  investigated  In  a more  conven- 
ient way;  in  particular,  the  general  theory  of  linear  Integral  equations  can 
now  be  brought  into  play  in  an  easier  fashion. 


3.  The  parametrix- 

We  have  to  study  the  character  of  the  singularity  of  g(P,Q; £ ) if  P moves 

& * 

to  Q.  We  refer  the  distance  function  r(P  ,Q  ) back  to  the  domain  B and  find 

(2.3.1)  r(P*(P),Q*(Q))  - (ST  (x*-  f*)2)372 

j-1  J ^ J 

L 3 *«,*<?,> 


cii  ‘tt(tj)(V  Si)(V  Tk’*2  af,  i,(xi-?i)(v 

♦0(r4»l/2  , 


r *»  r 


(psq) 


with 

(2.3.2) 


"5  * 

3x,  'dx. 


aik(xP  " ^ dx± 


By  virtue  of  the  continuous  differentiability  of  G(P  ,Q  ) - r(P  ,Q  ) 1 
and  the  character  of  the  transformation  (2„2),  we  can  assert  that 


(2.3.3)  g(P)Q?£  ) ■ — * — ~~  + h(P,Q;£  ) , 

4rrr(P  (P),Q  (Q)) 

where  h(P,Q;  £)  has  continuous  first  derivatives  in  D, 

1 * *\-l 

The  function  r(P  ,Q  ) has  for  P near  Q the  same  asymptotic  behavior 
as  g(P,Q| £ ) in  the  sense  that  the  difference  function  h(P,Q; £ ) is  continuously 
differentiable  in  D.  It  can  even  be  shown  that  the  second  derivatives  of 
h(P,Q;£)  become  infinite  at  most  like  r(PaQ)  ^ if  P approaches  Q.  We  can 


- 13  - 


construct  another  function  of  P and  Q which  is  somewhat  simpler  and  has  the 
same  asymptotic  behavior . In  fact,  the  function 


(2*3*4)  S£(P,Q) 


_1_  rl 

4tt  L2 


^aik^^  + aik^xj ^ ^xi~  ?k^ ' 


i,k-l 

satisfies  our  requirements,  as  is  Immediately  seen  by  series  development  near 
the  point  Q. 

We  define  now  the  parametrix  functions  s(P>Qs£,)  with  respect  to  the 
differential  equation  (2.2„9)  as  follows: 

a)  The  parametrix  s(P,Q;£)  is  a symmetric  function  of  both  arguments. 

b)  We  can  write 

(2.3.5)  s(P,Q;€  ) = Se(P,Q)  + R£(P,Q)  , 


where  R^(P,Q)  is  twice  continuously  differentiable  in  D for  P/Q,  has 
uniformly  bounded  first  derivatives  for  P^Q,  and  has  second  derivatives 
which  become  infinite  at  most  like  r(P»Q)  ^ for  P— ^Q. 

c)  The  parametrix  vanishes  if,  for  fixed  Q€  D,  the  argument  point  P 
lies  on  the  boundary  surface  C. 

The  construction  of  such  a parametrix  can  be  performed  in  various  ways. 
The  concept  of  the  parametrix  was  introduced  by  Hilbert  [11],  who  applied 
it  to  study  the  dependence  of  solutions  of  partial  differential  equations  on 
parameters  which  occur  in  the  coefficients  of  the  equation.  While  it  may  be 
difficult  to  obtain  a fundamental  solution  for  a given  differential  equation, 
the  construction  of  a parametrix  requires  no  comparable  efforts,  since  only 
the  boundary  condition  and  the  singular  character  have  to  be  observed. 


4-  Tlje  .lnt9g£al_§quati&as_fqx  g(P.sQv  £ )- 

Let  Q [u,v]  be  the  bilinear  form  belonging  tc  the  quadratic  functional 
(2.2.10).  Let  u and  v be  two  functions  which  are  twice  continuously 
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differentiable  in  D+C.  We  have  by  Green's  first  identity 

(2.4.1)  Q^[u,v]  * ,Vu)d<r  -|j|vLg[u]dr  , 


D 


where 


(2.4.2)  Pe(-^,Vu)  - XI  A,k 

and  "V  is  the  inner  normal  to  C with  components  From  this  result 

follows  Immediately  Green's  second  identity 

(2.4.3)  (vL^tul- ulgtvDdi:  a " »^u)  “ uP£(  V ,V  v))d<7" 

D C 

We  apply,  at  first,  (2„4.l)  to  the  two  functions 


(2.4.4)  u =*  g(P,Qj  6 ) , v » g(P,R?  SQ)  , Q,H€D 

We  have  to  modify  the  identity  in  the  usual  way  by  eliminating  from  the  domain 
of  integration  small  spheres  around  the  points  Q and  R where  the  integrand  is 
singular.  If  we  let  the  radii  of  the  two  safety  spheres  tend  to  zero,  the 
surface  integral  over  the  sphere  around  R will  disappear,  since  v tends  to 
infinity  only  like  the  reciprocal  of  its  radius,  while  its  surface  area  tends 
to  zero  like  the  square  of  its  radius.  The  contribution  of  the  sphere  around 
Q is  given  by 


(2.4.5)  g(Q,Rj  e0)  4rr 


Aik^  f,> 


dw 


U 


3/2 


i.k®! 


where  dw  is  the  surface  element  of  the  unit  sphere  U around  the  point  Q. 
From  definitions  (2.8)  and  (3.2)  follows 

3 

(2.4.6)  £ Vsj),kl(fj)  ' Si£°(V 


k-  1 
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There  is  no  restriction  of  generality  if  we  assume 
3 3 

(2.4,7)  XI  a (9  ) y±K  Ai  v\ 

i,k-l  1 1JC  i-l1  1 

that  is,  that  the  quadratic  form  is  at  Q on  principal  axes.  It  is  well  known 
from  the  theory  of  the  attraction  of  ellipsoids  that 


(2.4.8) 


JM  -I 


00 


[k  + k 1/*  + A -V ^ 

U L11  A2  2 3 3J 


an  \ 

J /V(A1+t)(A2+t)(A3+t) 


Hence,  we  obtain  easily 

(2-4.9)  t&.-l  ft 


J&i. 


V-l 


2 ^ [A  + A *2^  + A "l^/2 
U L 1 1 *2  2 A3  3J 


- «2TT‘(A1A2A;3) 


-1/2 


Thus,  we  have 

(2.4.10) 


<J-T 


J&L 


- e( 


i3/2 


-1  « 


and  the  contribution  of  the  sphere  around  Q becomes  simply  g(Q,R:  6 ) . 

0 

We  have  further  Lg[u]  =.0  in  D,  whence  we  derive  finally 


(2o4.ll) 


g(Q,Rj £ ) - Qf.[g(P,Qj^),g(P,Rj  £rt)] 


o'  6 

In  exactly  the  same  way,  we  find 

(2.4.11')  g(Q,R|  6)  » Qp  [g(P,Q;  £ ),g(P,R;  £ ) 1 

We  define  now  the  difference  terms 

(£,£  ) 


(2.4.12) 


A 


ik  ° 3 V^je)-Aik(x„j  ^ 


(£,  £ ) 


y ° - y^(x„;  £)-  /(x^  £Q) 

and  the  corresponding  bilinear  form 

(£,£  ) (£,£  ) ~ 

• ,5,  *» " n *k  • r 


(e<£o) 


uv 
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Subtracting  (4.1l)  from  U.ll'),  we  then  obtain 

)]d'C 

This  is  an  integro-differential  equation  for  the  Green’s  function  g(Q>Rj£  ) 
in  terms  of  the  Green's  function  g(Q,R;  E.  L It  is  very  simple  and  symmetric, 
but  has  the  disadvantage  that  the  coefficients  of  the  integro-differential 
equation  are  highly  singular  and  thus  this  equation  does  not  exhibit  clearly 
the  continuous  character  of  g in  dependence  on  S' . 

In  order  to  overcome  this  difficulty,  we  apply  now  (4.3)  with 

(2.4.15)  u»g(P,Q;£)  , v » s(P,Rj  £ ) , Q,RsD  , 

where  s Is  any  admissible  parametrix.  We  have  ..o  exclude  again  the  singular 

♦ 

points  Q and  R by  safety  spheres  and  to  pass  then  to  the  limit  of  vanishing 
radius.  As  before,  we  find  easily  that  the  contribution  of  the  spheres  around 
Cl  and  R are  s(Q,R;  £)  and  -g(R,Q;  £),  respectively.  Using  the  fact  that  g is 
a solution  of  (2.9)  and  is  symmetric  in  its  arguments,  we  thus  obtain 

(2.4.16)  g(Q,Rj  £)- s(Q,R?  6)  -|jJg(P,Qi  e)Le£s(P,Rj &)]drp 

D 

This  equation  is  an  integral  equation  of  the  second  kind  for  g(Q>Rj  C ) with 
an  admissible  kernel  for  the  Fredholm  theory.  In  fact,  we  may  put 

(2.4.17)  L8[s(P,R?£  )]  - L£[s(P,R;  £ )- g(P,R;  £.)]  , 

and  by  definition  of  the  parametrix  the  second  derivatives  of  s-g  become 
infinite  of  at  most  the  same  order  as  r(P,R)  if  P-^-R.  However,  we  want 
to  transform  the  integral  equation  into  such  a form  that  it  may  be  resolved  by 
a Neumann's  series. 

For  this  purpose,  we  apply  (4.3)  again  with 

(2.4.18)  u-g(P,Q?E) 


err  (e,ej 

(2.4.14)  g(Q,Rj  £ )-  g(Q,R;  e ) - - I E Cg(P,Q?  £ ),g(P»Rj  £ 


V - g(P,R;  £q)-  s(P,R;  £.q) 
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Observe  that  v has  now  bounded  first  derivatives  in  D,  even  near  the  point  R. 
Hence,  there  is  no  contribution  from  the  safety  sphere  around  R and  we  obtain, 
since  Lg[u] * 0 , 

(2.4.19)  g(Q,Rj  E^J-sCQjR;  €q)  * L£[g(P,E}  eo)-s(P,R;  £q) ]g(P,Q;  €.)d Cp  . 


Subtracting  (4.19)  from  (4. 16),  we  obtain 
(2„4»20)  g(Q,R;E)- g(Q,R;£o)  - s(Q,Rjfc)  - s(Q,R;Eo)  +, 


(e*e  ) 


(P,R)g(P,Q;e)drp  , 


with 

(2.4.21)  K 


(£,£•) 


(P,R)  - L.[g(P,RSe  )]-L  tg(P»R)E  )]♦  L-[s(P,Rjl)  - a(P,R;E  )] 


Up  to  this  point,  no  particular  assumption  was  made  concerning  the  dependence 
of  the  parametric  upon  the  parameter  £ . We  observe  that  the  a^Cx^jE)  depend 
analytically  upon  £ and  that  we  can  assume  without  loss  of  generality  that 
the  parametrix  depends  on  £ in  a sufficiently  smooth  manner.  Under  this 


assumption,  the  kernel  K 


is  small  of  the  order  £ - 6 , except  near  the 


point  R,  where  it  becomes  infinite.  But  we  may  write 

(£,£-) 

(2„4o22)  K ^ (P,R)  = Lp[g(P,R;£  )- s(PsR;£  )'J  + L-[s(P,Rj£)- g(P,R;£)]  , 

c.  o o o 

and  in  view  of  the  characteristic  property  of  the  parametrix  we  may  conclude 

(E  6 ) / ”1 

that  K ’ 0 (P,R)  becomes  infinite  at  most  like  r(P,R)  if  P~>R.  For 

|£ - £q|  small  enough,  the  Neumann's  series  for  the  reciprocal  kernel  of 

r(£,£0)  ... 

R u will  converge . 

Let  us  put 

(2.4.23)  $ (QjR;  £• , € ) * g(Q,R;  £■  ) + s(Q,Rj  £ ) - s(Q,R;  £q) 

We  have 


(2.4.24) 


(s.e  ) 


(Q,R)  - L„[»(Q,R;  E,  £ )1 


and  the  integral  equation 
(2.4.25) 


g(QsR;&)  » ^(Q,R5^£0)+JJl£[  ^(PSR5£,£  )]g(p,Q|Odrp 

D 

for  the  unknown  g(Q,R;£)„  This  integral  equation  is  of  the  Fredholm  type 
and  admits  resolution  by  a Neumann's  series. 
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We  might  have  chosen  s(P,Q$£)  to  depend  even  analytically  on  6 . This 
choice  is  of  interest  in  the  case  that  the  coefficient  p(x. ) of  (l.l)  depends 
analytically  on  its  variables.  In  fact,  we  can  th3n  assert  that  Kv  * ° 
depends  analytically  on  £ and  we  have  the  following  results  If  p(x^)  is  an 
analytic  function  of  the  x^'s,  the  Green3s  function  g(P,R;£)  depends  analytically 
upon  its  parameter  £ . 

Under  our  more  general  assumptions,  we  can  only  assert  that  g(Q,R; £ ) is 
continuously  differentiable  with  respect  to  £.  For  this  purpose,  we  have  to 
assume  only  that  the  parametrix  depends  different iably  on  £ . We  divide  the 
identity  (4-20)  by  ( £-  £^)  and  pass  to  the  limit  € ■ We  see  immediately 
that  the  partial  derivative  of  g(Q,R;  £ ) with  respect  to  the  parameter  £ exists 
and  satisfies  the  equation 


(2,4-26)  - ^ +jjjl>p[g(P,R5  £ ) ]g(P,Qj  £ )d*C 


at 


+IVlaafi]]g(p,Q5  £ )dr: 


where 


i,k*l  ^xi 


^ ~ 2A„. 

(2,4.27)  L'Cu]  - 


-2*3  - u 

ax,  J 


We  see  also  from  the  integral  equation  (4-20)  that  the  first  partial 
derivatives  of  g(Q>R; E ) depend  continuously  on  £ „ Using  this  fact,  we  can 
divide  the  identity  (4-14)  by  ( £ - £ ) and  pass  to  the  limit  £ 3 £ j we  find 


(2,4.28) 


3 £ 


E;[g(P,Q;  £ ),g(P,Rj  £)]dTT, 


where  Eglu,v]  is  the  bilinear  form 

(2  L 29)  E/[u  vl  - T -2-  A It  ■ t)  ^ -2i  . 

12.4.29;  E^LUjVj  Aik^,£;  ^ ££ 


i,k*i 


'dxs 


uv 
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Formula  (4-28)  Is  not  only  more  symmetric  than  (4.26),  but  it  hasy  moreover,  the 
advantage  of  being  independent  of  the  choice  of  the  parametrix.  It  is  a typical 
variational  formula,  expressing  the  change  of  the  function  with  the  parameter 
in  terms  of  the  function  itself  and  its  partial  derivatives  with  respect  to  its 
arguments . The  case  £ * 0 is  of  particular  interest  for  most  applications  and 
will  be  considered  in  detail  later  on. 


5 . Soas.  Jtogqa&3JL£iaa  • 

We  return  to  the  integral  equation  (4»20)  and  observe  that  the  resolution 
by  a Neumann’s  series  leads  to  a development  of  g(P,Rj£  ) in  terms  of  increasing 
order  in  (£-  £o).  We  obtain  thus  a numerically  convenient  formula  for 
determining  the  Green’s  functions  of  near-by  domains  with  an  arbitrary  degree 
of  approximation. 

We  want  to  point  out  one  feature  of  this  development  which  leads  to  interesting 
inequalities  in  an  important  special  case.  Iterating  (4.25),  we  obtain 


(2.5.1)  g(Q,R;£  ) - #(Q»R;  £. 


[tf  (P,R;  8,  E )]  y(F,Q;  £ , S M'T. 


Jl!Sr6tis'(p’R!E> 

D D 


£o)]g(P,Mj  £ )dTpdirM  . 


In  the  case  that  the  coefficient  p(x.)  of  the  differential  equation  (l.l)  is 

J 

non- negative,  we  observed  in  Section  1 that  G(P,M)  is  a positive -definite 
kernel;  by  definition  (2.3)  the  same  is  then  true  of  g(P,M;S).  Hence,  the 
last  double  integral  in  (5 . 1 ) is  non- negative  and  we  are  led  to  the  inequality 


(2„5.2)  g(Q,RjE  )>  ^(Q,Ri  £.,  £ ) +1 


L£[#(P,Rj  e,  eo)] 


(p, Q;  e. 


£o)d1:P 


We  can  generalize  this  last  inequality  as  follows.  Let  be  a set  of  N points 
in  D and  Xy  be  N arbitrary  real  numbers . From  the  positive-definite  character 
of  the  kernel  g(P,M;  6)  and  the  identity  (5«l),  we  can  then  derive  the  inequality 
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(2.5.3)  XT  [g(QiSQki£)- ^(Q.,Qk^  e,£o)]A1Xk^J]\L£Ehlhd'i:  , 

i n k^l  - . 


where 


N 

(2,5.3')  h(P)  » XI  y(F»Q45  S,  6 ) A 

i-  1 1 


These  inequalities  become  particularly  simple  if  we  choose  the  parametrix 
s(P,Qj  £.)  in  such  a way  that  s(P,Qj  £Q)  * g(P.tQj  gQ).  *n  fact,  in  this  case 
we  have  by  definition  (4 .>23) 


(2.5.4)  # (Q,R;  £,  Sq)  - s(Q,R;£) 

The  inequalities  (5.2)  and  (5.3)  are  sharp,since  we  obtain  equality  if  we 
choose  s (P,Qj  £ ) « g(P,Qj  £ ) . 

It  is  not  difficult  to  derive  these  inequalities  in  a direct  and  elementary 
way,  but  it  should  be  observed  that  they  are  only  one  set  of  a great  number  of 
estimates  provided  by  the  integral  equation  (4.20)  in  the  case  of  a positive- 
definite  Green's  function.  For  example,  we  can  iterate  (5.l)  again  and  can 
express  g(Q,Rj  £ ) in  terms  of  known  quantities  and  a remainder  integral 

(2.5.5)  ^|m^[^(P,R)]L£[^(M,Q)]L£[^  (T,P)]L£[j>f  (S,M)]g(T,S;  £ 

D D D D 

which  is  again  non-negative.  Iterating  the  integral  equation  for  g(Q,R$  £ ) 
in  this  way,  we  can  approximate  the  unknown  function  arbitrarily  and  with  a 
remainder  term  of  non-negative  value.  Each  approximation  leads  at  the  same 
time  to  an  inequality  for  the  desired  Green's  function. 


6 . Variational  formulas  and  variational  tensors . 

Thus  far  we  considered  the  parameter  £ in  the  transformation  (2.2)  as  a 
sufficiently  small  but  finite  quantity.  Re  shall  now  obtain  considerable 
simplifications  in  our  formulas  if  we  treat  £.  as  an  infinitesimal,  that  is, 
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if  we  retain  in  all  our  formulas  only  terms  which  are  of  the  first  order,  at 
most,  in  £ . The  formulas  thus  obtained  are  indeed  variational  formulas,  since 
they  determine  the  Green's  functions  of  domains  which  are  in  an  infinitesimal 
neighborhood  of  the  original  domain  D,  belonging  to  £ *0. 

ffe  calculate  at  first  the  coefficients  A^(x^;  g ) up  to  the  first 
order  in  £ . Since  * 

£>X. 

(2.6.1)  0 x " ^ij  + £ 


3S. 
1 


we  have  in  the  required  degree  of  precision 

3 


(2.6.2) 

and 

(2.6.3) 


@(x. ) * 1 + £ —l1  + o(  £ ) 


. » 0x„ 

Xs  1 X 


3S. 


c?x^ 


“ sij" £ 3i,  + o(e) 


3 


Hence,  using  definitions  (2.8)  and  (2.8!),  we  find 

3s. 


as, 


9s, 


(2.6.4)  Aik(x.;  £ ) - ^ik+  S / ^ik  ^ ~ 3% 


~^l  + o( 


e) 


(2 „6. 4' ) J)(x.;  £ ) « p(x^)  + £ ^ I (pS^)  + o(  £ ) 


3 


J J 3*1’  ”J 

Now  we  are  ready  to  derive  an  elegant  variational  formula  from  the 
identity  (4.14)*  We  define  the  "variational  tensor" 


(2.6.5)  Tik(P;Q,R)  - 


+ ?.g(LtSjLal 


“ V.g(P?Q;o)  * Vpg(F,R;o)  S 


ik 


and  using  xt,  we  put  (4.14)  in  the  form 

\\\\^  3s  (p) 

(2.6.6)  g(Q,Rj  &)»  g(Q,Rjo)  - E \.\H  4—  ^(PjQ,*) 

“le 


I i.:k“l 


° YZ  (psv)g(p>Q5c)g(P,Ri  .)7dT0+  o(C  ) 

k-1  J 
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The  variational  tensor  T^(P;Q,R)  will  play  a central  role  in  the  sequel. 
It  is  based  only  on  the  original  Green's  function  g(P,Q;o)  and  is  to  be 
considered  as  known.  The  tensor  is  symmetric  in  its  indices  as  well  as  in 
its  dependence  on  Q and  R.  For  Q-R  the  tensor  becomes  a Maxwell  tensor  of 
the  type  frequently  used  in  electrostatics. 

The  trace  of  the  tensor  is 

(2.6.7)  T(PjQ,R)  = H T (PjQ,R)  - Vpg(P,Q5o)- Vpg(P,R;o) 

k- 1 r 

ffe  verify  easily  the  following  identity s 
3 

(2.6.8)  YZ  T,.  (P;Q,R)T.,  (P;Q,R)  - 

k-1  1 J 


^tTCPjQ.QjTy  (PjR,R)  + T(PjR,B)T  (PjQ,Q)  ] 


TCPjQjR^^PjQjR)  + $ T(PlQfQ)T(P,BfR) 


which  reduces  for  Q - R to 
(2.6.9) 


k-1 


Tik(F»«.<i)TjkCl’S<3,Q)  - SjjTtPsU.Q)2 


This  shows  the  important  fact  that  the  rows  and  columns  of  the  Maxwell  tensor 
T.^(P|Q,Q)  are  orthogonal  to  each  other. 

So  far  we  used  only  the  formal  structure  of  the  variational  tensorj  now 
we  utilize  the  fact  that  the  function  g(PsQjo)  satisfies  the  differential 
equation  (l.l)  in  dependence  on  each  variable.  We  find  then  easily 


(2.6.10) 


’ST"  tg(P,Qjo)g(P,Rjo)] 
k-1  J *xi 


This  simple  differential  identity  permits  us  to  bring  the  variational 
formula  (6.6)  into  t,he  simple  form 
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(2.6.11)  g(Q,Rj Z ) - g(QjR;o)  » 

ej|fl  <ZJ  tllk-.  p<p)«<F,Q„>g(p,R,o>  S^jtfllar^ote) 

We  have  thus  expressed  the  first  variation  of  g(Q,R)  as  an  integral  over 
the  domain  D whose  integrand  is  a divergence  term.  We  assumed  D to  have  a 
smooth  boundary  surface  C and  may,  therefore,  transform  this  integral  into  a 
surface  integral  over  C.  We  make  use  of  the  fact  that  on  C 

(2.6.12)  g(PyQ;o)  - 0 , Vpg(P,Q5o)  - ^ , P€.C  . 


We  observe  further  that  the  components  of  the  tensor  T^(PjQ,R)  become 
infinite  for  P“Q  and  P»R.  Hence,  we  have  to  exclude  these  singularities 
from  the  domain  of  integration  by  infinitesimal  spheres  and  to  take  into 
account  their  contributions.  After  an  easy  calculation,  we  find 

(2.6.13)  g(Q,R}£)- g(Q,R;o)  - -J  2-dZJL\<il  (£S*^)d <7* 

C P P 

3 

+ - ZI  SV(Q)  * sv(R)>  o(  £)  , 


k * 1 


2r, 


where  S is  the  vector  field  with  components  and  where  the  q^’s  and  r^'s 
are  the  coordinates  of  Q and  R,  respectively. 

It  becomes  now  convenient  to  return  to  the  original  Green's  function  by 
means  of  the  correspondence  (2.3).  We  denote  by 

(2.6.14)  £.  =*  S 7/ 

the  normal  shift  of  C under  the  infinitesimal  deformation  (2.2)  and  obtain 
up  to  higher  order  terms  in  & 


(2.6.15) 


S g(q,r) 


TjjGjPi&l 
JJ  92/v 

C 


2_G.QL.Rl 

-dv* 


S Z'pdCTp 
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This  elegant  variational  formula  was  derived  by  Hadamard  [10  ] in  the  case  of 
Laplace’s  equation  for  any  number  of  independent  variables.  The  more  complicated 
formulas  (6.6)  and  (6.1l)  have  the  advantage,  however,  of  being  valid  also  for 
domains  whose  boundary  surface  C is  no  longer  smooth.  We  shall  discuss  this 
extension  of  the  formulas  in  the  next  section.  Hadamard' s formula  is  frequently 
used  in  variational  considerations  because  of  its  great  formal  simplicity ; it 
leads  often  to  a heuristic  solution  of  extremum  problems  which  must  then  be 
established  precisely  by  a finer  technique  which  can,  in  general,  be  based  on 
the  formulas  of  interior  variation  of  the  type  (6.6). 

We  have  connected  in  this  ->®ction  Hadamard’s  variational  formula  with  the 
theory  of  the  variation  of  the  Green's  function  in  a fixed  domain  under  the 
variation  of  a parameter  in  the  coefficients  of  the  equation.  This  latter 
theory  i3  essentially  due  to  Hilbert;  it  provides  a simple  proof  for  the 
Hadamard  formula  and  permits  an  evaluation  of  the  error  term  arising  from 
neglect  of  the  higher  order  terms. 

7 . garitenatan  the-gagiatlQ.hal  formula  • 

It  is  obviously  necessary  to  extend  the  variational  formula  (6.6)  to  the 
most  general  domain  C in  space  for  which  a Green’s  function  exists.  In  fact, 
if  we  are  dealing  with  extremum  problems  for  domains  D involving  their  Green's 
functions,  we  will  have  to  characterize  the  extremum  domain  by  varying  it  and 
comparing  its  Green's  function  with  that  of  near-by  domains.  In  this  way,  we 
will  be  able  to  express  analytically  its  extremum  property  in  the  form  of 
identities.  But  we  are  not  sure,  a priori,  that  the  boundary  surface  C of  the 
extremum  domain  D is  smooth;  hence,  we  cannot  apply  the  results  of  the  preceding 
section  without  getting  rid  of  the  assumption  of  smoothness  on  C.  We  will  show 
now  in  this  section  that  the  variational  formula  (4-14)  holds  in  the  most 
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general  case  and  that  the  formula  (6.6)  which  has  been  derived  from  it  in  a 
formal  way  is,  consequently,  always  applicable. 

Let  D be  an  arbitrary  domain  in  space  which  possesses  a Green’s  function 
G(P,Q)  with  respect  to  the  differential  equation  (l.l).  It  is  easy  to  see  that 
the  deformations  (2,2)  carry  D into  a new  domain  D which  will  also  have  a 

Green's  function  G (P,Q)  if  £ is  small  enough.  We  can  then  define  in  D the 

functions  g(P»Q;£.)  as  we  did  in  Section  2. 

Let  D^  be  a sequence  of  domains  with  smooth  boundary  surfaces  which 

converge  to  D in  such  a way  that  D^CD^CD,  We  denote  by  gn(F,Q;  £ ) the 

Green's  function  of  corresponding  to  g(P,Q;£).  We  can  define  gn(P»Qj£) 
as  a piecewise  smooth  function  in  D by  putting  g^(P,Qj  £)*0  if  either 
argument  point  lies  in  D-D  . Putting 

(2*7.1)  u » g(P?Q;  £)  V « g^CPjRj  £q)  , Q,RSDn  , 

we  can  now  obviously  apply  the  first  Green's  identity  (4.l),  and  taking  notice 
of  the  singularities  of  u and  v,  we  find 

(2.7.2)  g^(Q?R;  6q)  * QgCg(P>Q;  £ ) ,gn(PjR;  £o)] 

Now,  we  can  pass  to  the  limit  n-— >oo;  because  of  the  well-known  continuity 
property  of  the  Green's  function,  we  obtain 

(2.7.3)  g(QsR;  £ ) » QFtg(P,Qj  £ )»g(F,R;  £ )] 

This  is  formula  (4.1l),  valid  now  for  the  most  general  domains  which  have  a 
Green's  function  at  all.  In  the  same  way,  we  can  prove  the  analogous  identity 
(4.11 1 ) and  the  desired  formula  (4»14)  follows  again  by  subtraction  of  the 
two  identities.  Thus,  the  variational  formula  (6.6)  has  been  established  in 
the  most  general  ease. 

It  is  to  be  observed  that  we  are  using  here  the  terra  Green's  function  in 
the  generalized  sense,  namely,  as  the  limit  of  the  Green's  functions  (in  the 
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proper  sense)  of  smoothly  bounded  exhausting  subdomains.  In  general,  it 
cannot  be  asserted  that  the  generalized  Green's  function  vanishes  at  all 
points  of  the  boundary  G of  D.  ®ut  it  is  necessary  to  use  this  concept  if 
we  want  to  establish  a useful  functional  analysis  with  respect  to  the 
differential  equation  (l.l). 

Our  preceding  results  are  in  no  way  restricted  to  the  case  of  three 

independent  variables,  but  can  be  generalized  to  n variables  if  the  nature  of 

the  parametrix  singularity  is  properly  adapted.  The  case  n»2  is  of 

particular  interest,  since  the  use  of  complex  variables  permits  various 

interesting  simplifications  in  the  formulas . 

We  start  again  with  the  differential  equation 

2 2 

(2.8.1)  LCu]  * V2u~  pu  - 0 , V2  - XT  -2-r 

i-l  9*1 

and  we  call  G(P,Q)  the  Green’s  function  of  (8.l)  for  a plane  domain  D if  it 
satisfies  the  three  requirements 

a)  For  fixed  Q£D,  G(P,Q)  is  a solution  of  (8.l)  as  a function  of  P. 

b)  G(P,Q)  +^f  log  r(P,Q)  is  continuous  in  D,  has  uniformly  bounded 

first  derivatives,  except  possibly  at  Q,  and  its  second  derivatives  grow  at 
_ ’1 

most  like  r(P,Q)  if  P approaches  Q. 

c)  G(psQ)  vanishes  for  P£C,  Q6D,  where  C is  the  boundary  curve  of  D. 
We  introduce  again  infinitesimal  deformations  based  on  a vector  field 

S (x.)  (i,j*ls2)  and  we  introduce  Green's  functions  g(P,Q;E)  in  the  same 

way  as  in  Section  2.  We  obtain  again  the  variational  formula  (6.6)  with  a 
variational  tensor  (6.5);  the  only  difference  is  that  the  indices  i,k  range 
only  over  1,2  instead  of  1,2,3. 
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There  occur,  however,  as  if  by  accident,  remarkable  relations  in  this 


special  case,  We  find 
(2.8.2)  TU(P;Q,R)  = 


^x. 


l°2  _2si£jiLial  _ /?g(PaQ;o)  ffgCP.tRjol 


9x„ 


T22(P;Q,r) 


t (p.Q  r)  • 3 g(P*ftsft)-  ,9-gif^-iol  + gg(Pj.Q;ol  ^(P.Rjo). 

i12'r,^s'tt^  ^x_  ^x_ 


3x, 


■25  x. 


T21(P;QsR) 


Thus,  T^  is  a symmetric  orthogonal  tensor  with  trace  zero;  we  have 
2 

(2.3.3)  5ZI  TikTJk  “ S±i  lvpg(P,Q;0)|2lVpg(P,R;o)|2 


We  introduce  the  complex  differential  operators 

d 

- 1 , — - - v— 

l 

and  we  write  (with  PSz) 


(2.8.4) 


-2.  .1  i-3- 

7$  z 2 ' 2x, 


-2_  - I (-iL_  + i -2-) 

' 'O-  1 ' > 


3x  ' * 3“  2 v 2x.  * 0X.. 

2 £?  Z 1 << 


Xi+iX2  ’ 


(2.8.5) 


TU(P;Q,R) 


T12(P;Q,R) 


Re- 


&eS2j&is1  ^gCP^Rto) 

?z  3 Z 

lm  C?.k( M;o)  3g(P.R;o) 

| 3 Z ^ 25 


- T22(P;Q,R) 


56 


For  the  sake  of  uniformity,  we  collect  also  the  deformation  vector  field 
S^Xj)  into  one  complex  function 

(2„8.6)  F(z,z)  * S^(x^)+i  S2(x^) 

The  points  P,Q,R  are  replaced  by  their  complex  coordinates,  z,^, 
respectively.  With  these  notations, we  can  cast  (6.6)  into  the  complex  form 

(2.8.7)  g(^,  rjj£)- g(?,  *]5o)  » £ Re f j|[ 3 ^ ^ 3lg ^ 

UJ 

- 2g(z,  ^ ;o)g(z,  ^ ;o)  ^ 


+ o(  £ ) 
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ffe  recognize  that  the  variational  formula  becomes  particularly  simple  if 

F»F{z)  is  an  analytic  function  of  the  complex  variable  z.  In  this  case 

- 0 and  (8.7)  reduce.s  to 

^z 

(2.8.8)  g(  'Jj  £)- g(^,  ^;o)  - -2£  jJg(z,^;o)g(z,  ^;o)Ee^-^  (pF)J-d"Cz 

+ o(&)  . 

If  we  assume,  in  particular,  p =;0,  we  find  8g(^,  “ 0>  which  is  the  well- 

known  invariance  of  the  Green's  function  of  Laplace's  equation  under  a conformal 


mapping . 

We  can  obtain  from  (8.7)  a variational  formula  for  the  original  Green's 
function  G(p,q)  connected  with  the  fixed  differential  equation  (8.l)  and  the 
varying  domains  D . Let  G(P,Q)  denote  the  Green's  function  of  the  original 
domain  and  let 


(2.8.9)  8 G(P,Q)  - G*(PsQ)-G(P,Q)  + o(8  ) 

denote  the  first  variation  of  this  Green's  function,  if  the  domain  D is 
transformed  into  D by  a deformation  (2.2).  Then,  we  find  easily  from 
(8.7)  and  Taylor’s  theorem 


(2.8.10) 


£8(9,  r) 


> - £Be 


s 


[8 


'B  z 


B z 


D 


2z 


- 2G(: 


9 . Singular  _variatioixs_lnJih^-Plane  • 

We  saw  in  the  preceding  section  that  interior  variations  based  on  an 
analytic  function  F(z)  lead  to  a particularly  simple  variational  formula  for 
the  Green's  function.  In  the  very  important  special  case  of  the  Laplace 
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equation  in  the  plane,  which  is  closely  related  to  the  theory  of  analytic  functions, 
the  simplification  becomes  even  too  great j since  the  Green’s  function  is 
invariant  under  conformal  mapping  in  this  case,  we  find  only  that  the  first 
variation  of  g(P,Q;£)  is  zero.  In  order  to  obtain  a more  applicable  variational 
theory  of  conformal  mappings  and  univalent  functions,  it  proves  useful  to 
consider  a slightly  more  general  class  of  deformations.  These  deformations 
are  based  on  a complex  function  F(s,z)  which  is  analytic  in  the  domain  D 
considered,  except  in  a small  fixed  circle  inside  of  D,  The  same  type  of 
variation  is  also  of  interest  in  the  case  of  more  general  differential 
equations,  as  will  be  seen  from  the  applications  in  Chapter  IV. 

Let  zq  be  a fi’ed  point  in  D and  let  K(E,  ^z^)  denote  the  circle 
I z-  zq]  V E « p.  We  assume  p so  small  that  Kc.D  and  define 


Ijbl 


i<* 


(2.9.1)  F(z,z) 


6 z - z 


, i 

J-2—  (6 

6 z - z 


i Z-2 


for  zeD-  K 


+ 3 


z-z  ] 
* 


•) 


for  z£K  . 


V_  “ P P jr 

By  easy  calculations  it  can  be  seen  that  this  deformation  field  is  twice 

continuously  differentiable  in  Dj  it  yields  a one-to-one  mapping  of  D into 

a schlicht  near-by  domain.  Thus,  our  general  theory  of  variations  applies. 

We  insert  (9 . l)  into  (8.10)  and  find  easily 

(2.9.2)  e1*  ~ 

-2P(^,^)  ie1Kjtt(Z(r)0(z,y)  -2.  (J^)atr 


Let  us  suppose  that  the  argument  points  ^ and  to  lie  outside  of  the 
circle  K(£.,zo).  Formula  (9.2)  represents  then  the  first  order  change  of  the 
Green’s  function  with  respect  to  D and  for  the  differential  equation  (8.1), 
if  the  boundary  C of  D is  mapped  by  means  of  the  analytic  function 
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(2.9.3) 


, i<*  .2 

f (z)  * 2*7  -ft-  ' 
6 z “*  z 

o 


which  is  univalent  on  C,  this  mapping,  G is  transformed  into  a new  curve  C 
which  bounds  a new  domain  D and  (9.2)  gives  the  change  of  the  Green’s  function 
under  this  variation. 

Variations  of  the  type 


(2.9.3') 


z + 


-P- 


Z-  7, 


have  been  used  successfully  in  the  theory  of  schlicht  functions  and  are 
called  "interior  variations"  of  the  domain.  The  factor  ^ which  occurs  in 
(9.3)  is  unessential  and  has  only  been  introduced  in  order  to  permit  a simple 
continuation  of  the  deformation  into  the  critical  circle  K(£.,zq)  which 

remains  one-to-one  inside  this  circle.  We  reformulate  our  result  (9.2)  in 

. . * 
the  following  ways  Let  (9.3 5 ) transform  the  domain  D into  & and  map  the 

points  ^ and  rj  into  ^ and  fj  . If  G (^,  1^)  is  the  Green’s  function  of  the 

new  domain  D , we  have 

1 ott  o_J> aii__ 


(2.9.4)  G (|  , (J  )-G(^,  rj)  - Re|8TTjCy 


- 2 p2eicK 


JJ  0lz- 


9 z 9z 

o o 


D 


This  result  is  obviously  of  particular  elegance  in  the  case  of  the  Laplace 
equation  (p=0),  but  is  also  useful  in  the  more  general  case.  Its  main 
advantage  lies  in  the  fact  that  no  derivatives  of  the  Green’s  function  occur 
under  the  integral  sign. 


10. 


Neumann's  function . 


We  derived  in  Section  6 an  interior  variational  formula  for  the  Green’s 
function  by  using  the  identity  (4.14)  and  the  fact  that  g(P>Qj£)  and  its  first 
partial  derivatives  depend  continuously  on  £ » We  can  derive  in  the  same  way 
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analogous  formulas  for  various  other  important  domain  functions  related  to  the 
differential  equation  (l.l). 

We  define  the  Neumann's  function  N(P,Q)  for  a three-dimensional  domain  D 
with  respect  to  the  differential  equation  (l.l)  as  usual* 

a)  N(P,Q)  is,  for  fixed  Q£.D,  a solution  of  (l.l)  as  a function  of  P. 

b)  It  has  In  D and  at  Q the  same  differentiability  properties  as  the 
Green's  function. 

c)  For  fixed  QeD  and  P^C,  we  have 


(2.10.1) 


ggttt .ft)  . o 


03/ 


Let  us  assume  that  for  a given  domain  D such  a Neumann's  function  does 
exist.  We  introduce  the  infinitesimal  deformations  (2.2)  and  obtain  near-by 
domains  D with  Neumann's  functions  N (P,Q),  if  )£|  is  small  enough.  Let  us 
define  in  analogy  to  (2.3)  the  functions 

(2.10.2)  n(P,Q;e)  - N*[P*(P),Q*(Q)] 


in  the  fixed  domain  D,  and  let  us  study  their  dependence  on  6 . We  observe 
that  n(P,Q; 6 ) is  symmetric  in  P and  Q,  in  view  of  the  well-known  symmetry 
of  the  Neumann's  function;  it  satisfies  in  each  variable  the  differential 
equation  (2.9),  since  N (P,Q)  satisfies  (l.l)  in  D . 

In  order  to  establish  the  conditions  which  n(P,Q; £ ) fulfills  on  the 
boundary  surface  C,  we  make  the  following  observation.  Let  U (P  ) be  an 

arbitrary  solution  of  ^1,1)  in  D which  is  continuously  differentiable  in 

* * * 

D +C  ; let  V (P  ) be  continuously  differentiable  in  L +G  , We  define  in  D 

the  functions 

(2.10.3)  u(P)  » U*[P*(P)1  , v(P)  » V*[P*(P)] 


and  find  by  (2.6)  and  (2.10) 
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(2.10.4)  jJJ  Q£[u,vld'C  - jjj  [V  U*  • V V*  + pU* V* ]d^* 

D D* 

We  use  now  Green’s  first  identity  on  both  sides  of  (lQ.4)  and  observe  that 

-}£  "N* 

u satisfies  in  D the  equation  (2.9),  while  U satisfies  in  D equation  (l.l). 
Thus,  we  obtain 

(2.10.5)  ^ vP£t>,Vu]d<r  = JSV*  -^hrd <T* 

C c*  dv 

Since  v(P)  is  arbitrary,  we  obtain  from  this  relation  the  identity 

* 

(2.10.6)  Pc[^,Vuld(T  » dCT* 

which  clarifies  the  meaning  of  the  important  linear  functional  P£ . We 
obtain,  in  particular,  the  following  boundary  condition  for  the  function 
n(P,Q;£  ): 


(2.10.7)  PgCv,  ^7pn(P,Q;  Z )]  - 0 , for  P6  C,  Q€D  . 

If  we  now  apply  (4°l)  with  respect  to  two  functions 

(2.10.8)  u * n(P,Q;  C ) , v = n(P,R;  E)  , Q,R£D  ; 

we  obtain  by  the  same  calculations  as  were  already  performed  in  Section  4 


(2.10.9)  n(Q,Rj  £q)  » Qg(n(P,Qj  £ ),n(P,Rj  £q)] 

Interchanging  t with  and  subtracting  the  resulting  formula  from  (10.9), 
we  prove  finally  the  identity 


(2.10.10)  n(Q,R; £ ) - n(Q,R;  £ ) 

o 


(E,€ . ) 

E 0 Cn(P,Qj  £ ) ,n(P,Rj  £ ) ]<fCp 


We  do  not  enter  here  into  the  detailed  proof  that  n(Q,R;£  ) and  its  first 
partial  derivatives  depend  continuously  upon  the  parameter  £ . We  proceed 
immediately  to  derive  from  (10.10)  in  a formal  way  the  variational  formula 
for  N(Q,R) * n(Q,R;o) . We  define  the  variational  tensor  for  the  Neumann’s 


function  to  be 
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(2.10.11)  r jpjq.b)  - -2-“  amM * m&al  a^Ri 

IJc  9 3 xi 

- VpNdSQj-VpNCP ,R)^ik  , 

and  we  derive  from  (10 .10)  the  result 

,~r  „ 

(?) 


(2.10.12) 


n(Q,R;£  )-N(Q 

3 


•*>  • #fe 


- XI  (pSk)N(P,Q)N(P,R)j  d^  + o(  £ ) . 

This  is  the  fundamental  interior  variational  formula  for  the  Neumann's  function, 
which  stands  in  complete  analogy  to  formula  (6.6). 

The  tensor  T^  satisfies  the  same  differential  relation  as  the  variational 
tensor  T^>  namely 

(2.10.13)  ^ ~£r  (P?Q,R)  - P(P)  (»(P,Q)N(P,R))  ; 

k-i  axk  ik  ?xi 


thus,  we  may  put  (10.12)  into  the  form 

(2.10.34)  n(Q,R;  Z ) - N(Q,R)  - ejf  (XI  &ik~  p(?)N(P}Q)N(P,R)  8^ 


_3  ± 

!W  axk  'i^l 


*si(p))cfrp  + o(e.)  . 

Following  through  the  same  calculations  as  in  Section  6,  we  obtain,  by 

&N 

integration  and  use  of  the  boundary  condition  “^7  * 0 on  C,  the  final  formula 

(2.10.15)  SN(Q,R)  iVpNCPjQj^Vp^PjRj  + ptPjHtP^jNfP^JlSrpdCTp  . 

C 

This  formula  corresponds  to  the  Hadamard  variational  formula  (6.15)  for  the 
Green's  function  and  plays  a similar  role. 

We  want  to  remark  finally  that  the  reasoning  of  this  section  can  easily 
be  extended  to  the  so-called  Robin's  functions  of  a domain  D with  respect  to 
the  differential  equation  (l.l).  A Robin's  function  is  defined  like  the 


I — 1 
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Neumann's  function  and  the  Green's  function;  but  instead  of  the  boundary 


conditions  G » 0 or 


3N 

2V 


0,  we  have  for  a given  positive  function  on  C 


(2.10.16)  - p-(P)R(P,Q) 


for  P€  C,  QeD  . 


In  order  to  carry  out  the  preceding  reasoning  in  the  simplest  way, 
it  is  convenient  to  vary  the  coefficient  |A(P)  together  with  the  domain  D 
in  such  a way  that  the  function 

(2.10.17)  r(P,Q;E)  =■  R*[P*(P),Q*(Q)] 


satisfies  on  the  boundary  surface  C the  condition 

(2.10.18)  £.)]  - |4P)r(P,Q;£) 

In  this  case,  it  is  easily  seen  that  the  identity  (4.14)  holds  also  fop 
r(P,Q;  E ) . Hence  we  can  obtain  the  same  formalism  as  before.  On  the  other 
hand,  it  is  not  difficult  to  derive  variational  formulas  for  the  Robin's  function 
of  a fixed  domain  and  a fixed  differential  equation,  but  with  varying  coefficient 
|jl(P).  In  this  way,  the  most  general  variational  result  may  be  obtained. 

11 . The  harmonic  Neumann's  function . 

The  results  of  the  preceding  section  are  not  immediately  applicable  to 
the  case  of  the  Laplace  equation 
(2.11.1)  V2u  - 0 , 


since  in  this  case  no  finite  domain  possesses  a Neumann's  function  in  the 
above  sense. 

It  is  customary  to  define  the  Neumann's  function  N(P,Q)  in  the  harmonic 
case  by  the  following  four  requirements! 

a)  N(P,Q)  is  harmonic  for  Pe  D,  except  for  a pole  at  Q. 

b)  N(F,Q)  - T3*-  r(F,Q)  1 is  harmonic  at  Q. 

4TT 
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c)  For  P€C,  Q6D,  N(P,Q)  has  a constant  normal  derivative, 

d)  We  have  the  relation 

(2.11.2) 

C 

The  last  condition  is  necessary  if  v?e  want  to  ensure  the  symmetry  of  the 
Neumann’s  function  in  dependence  on  its  two  argument  points.  We  have  by 
Green’s  theorem  the  identity 


5 


N(P,Q)dcrT3  =0 


(2.11.3) 


gJKEtfll 

2>v  P 


d<7"p 


9 


and,  since  the  Neumann's  function  has  a constant  normal  derivati/e  for  fixed  Q, 
we  can  sharpen  condition  c)  to  the  statement 

(2.11.4)  “ 4(C)  ^ , 4(C)  » surface  area  of  C. 

This  result  shows  that  the  constant  value  of  the  normal  derivative  is 
independent  of  Q and  that  for  any  two  points  Q and  R 


(2.11.5)  “^T  CN(P,Q)-N(P,R)]  - 0 , P€C,  Q,R€D  . 

The  function  N(PjQ,R)  = N(P,Q)  - N(P,R)  is  the  nearest  approach  to  the  general 
Neumann’s  function  concept  of  the  last  section.  It  has  only  the  inconvenience 
that  we  need  for  its  definition  the  normalization  (11.2)  for  its  components, 
which  cannot  be  extended  to  domains  with  more  complicated  boundary  surfaces. 

We  introduce,  therefore,  the  function 


(2.11.6)  N(P,WjQ,V)  » N(P;Q,V)- N(W;Q,V) 

* N(P,Q)  + N(W,V)  - N(P,V)  - N(W,Q) 

* N(Q?V;P,W) 

which  is  symmetric  in  argument  and  parameter  points  and  satisfies  the  following 


four  conditions J 
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a)  N(P,W;Q,V)  is  harmonic  in  P and  W,  except  at  Q and  V. 

b)  N(P,W;Q,V)- [r(P,Q)”1+r(W?vr1- rCPjV)"1- r(W,Q)-1]  is  regular 
and  harmonic  at  Q and  V, 

c)  N(P,W;Q,V)  has  for  fixed  Q,VeD  the  normal  derivative  zero  on  the 
boundary  surface  C. 

d)  H(P,P;Q,V)S0  . 

It  is  easily  seen  that  these  four  conditions  determine  N(P,W;Q,V)  in  a unique 
way  and  may  be  taken  as  its  definition.  Most  formulas  of  potential  theory 
become  more  elegant  if  this  function  is  used.  In  the  case  of  two  independent 
variables,  the  corresponding  expression  is  closely  related  to  the  symmetric 
integral  of  the  third  kind  on  a closed  Rieraann  surface,  obtained  by  completing 
the  domain  D considered  through  addition  of  its  double. 

We  might  obviously  define  a function  N(P,W;Q,V)  by 

(2.11.6' ) N(P,WjQ,V)  = N(P,Q)  + N(W,V)  - N(P,V)  - N(ff,Q) 

also  in  the  case  of  the  general  differential  equation  (l.l),  although  its 

role  is  rather  unimportant  in  this  case.  Let 

(2.11.7)  n(P,W;Q,V; £ ) - N* (P*,W* jQ* ,V* ) 

. * 
be  the  corresponding  solution  of  (2,9)  obtained  by  referring  the  function  W 

of  the  varied  domain  D back  to  D.  We  introduce  the  function 

(2.11.8)  N(P;Q,R)  * N(P,Q)  - N(P,R) 


and  the  variational  tensor 


ru  . ^ . . rv*  V-/  /-vJ 

Tik(P;R,W;Q,V)  = Tik(P;R,Q)+Tik(p;W,V)-Tik(p;R,V)-Tik(PjW>Q) 


N(P5R,W)  -f-  N(P;QSV)+-^-  N(P;R,W)  N(P;Q,V) 

axi  \ '*k  ^xi 


VpN(P°R,W).  ^pN(P;Q,V)  <5. 


ik 


(2.11.9) 
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We  then  deduce  immediately  from  (10.12)  the  re3U.lt 
(2.11.10)  n(R,W;Q,V;  £)- N(R,W;Q,V)  - 


In  this  formula , we  may  pass  to  the  limit  p = 0 because  of  the  existence 
of  a Neumann's  function  N(R,W$Q,V)  in  the  case  of  the  Laplace  equation.  We 
then  obtain 


(2.11.11)  n(R,W;Q,V;£  )-  N(R,WjQ,V)  - 

This  formula  is  the  basic  interior  variational  formula  for  the  Neumann's 
function  of  Laplace's  equation.  It  may  be  extended  to  the  most  general 
domains  D for  which  the  usual  generalized  Neumann's  function  exists. 

We  derived  in  this  section  the  variational  formula  (ll.ll)  by  a limiting 
process  from  the  corresponding  formula  in  the  general  case  (l.l).  It  is 
obvious  that  this  procedure  was  only  followed  for  the  sake  of  brevity;  one 
may  derive  (ll.ll)  directly  by  the  preceding  method  without  such  limit 
considerations , 

We  may  derive  from  (ll.ll)  and  the  equation 

(2.11.12)  ^ T (P;R,WjQ,V)  - 0 

k=  1 ?xk  lJc 

the  Hadar.iard  type  variational  formula 

(2.11.13)  ^N(R,W;Q,V)  - VpN(P;R,W)  • VpN(P;Q,V)^pd(7-p  . 

G 

It  should  be  observed  that  for  arbitrary  T£-D,  we  have 


(2.11.14)  V^(P,TjR,W)  - VpN(PjR,W)  , 

so  that  the  formula  (11.13)  might  as  well  be  formulated  in  terms  of  the 
symmetric  Neumann's  function  N(F,T;R,W)  alone. 
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12. 


We  shall  now  apply  the  results  of  the  preceding  sections  in  order  to  derive 
variational  formulas  for  some  quantities  of  physical  significance.  We  consider 
the  case  of  the  Laplace  equation  and  put 

(2.12.1)  G(P;Q)  » " g(P,Q)1 

Since  - G(P,Q)  is  regular  harmonic  in  the  domain  D and  positive  on  its 

boundary  surface  C,  we  have  by  the  maximum  principle 

(2.12.2)  g(P,Q)^0 

The  quantity 

(2.12.3)  C(Q)  - g(Q,Q) 

is  called  the  capacity  constant  of  the  domain  D (or  it3  boundary  surface  C) 
with  respect  to  the  point  Q.  Its  negative ,-C (Q) , is  the  potential  at  Q of  the 
charge  distribution  which  is  induced  on  the  conducting  surface  C by  a unit 
charge  placed  at  the  point  Q„  We  obtain  from  Hadamard's  formula  (6.15) 

(2.12.4)  5C(Q)  = 4TT  \ J (•a^^)2  Sp?d(r?  . 

C P 

This  formula  shows  the  monotonic  dependence  of  the  capacity  constant  upon  the 
domain  and  can  be  used  to  obtain  a clear  insight  into  the  dependence  of  G (Q) 
on  the  boundary  C.  Of  course,  one  can  find  also  an  expression  for  SC(Q)  in 
terras  of  a domain  integral  involving  the  tensor  T . Since  this  formula  is 
not  very  simple,  we  restrict  ourselves  to  one  important  special  case. 

We  suppose  that  the  domain  D is  the  exterior  of  a surface  G.  The  formulas 
of  Csction  6 are  also  valid  in  this  case.  If  Q lies  near  infinity,  we  may 
develop  the  Green' 3 function  in  terms  of  the  coordinates  q^  of  Q to  obtain 

. Jz,  2 7 

’ /’“VV 


(2.12.5) 


G(P,Q) 


* ort) 


P 


q3 
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Here,  G(P)  is  a harmonic  function  of  P in  D except  at  infinity}  it  vanishes  for 
P £ C.  It  is  easily  seen  that  G(P)  has  near  infinity  a development 

(2.12.6)  G(F)  - ^ El- * 0(^)1 

C(oo)  is  called  the  capacity  of  the  surface  C.  We  may  interpret  G(P)  as  the 
Green's  function  for  a point  charge  at  infinity}  this  notion  can  also  be 
justified  by  considering  the  transformation  of  G(P)  under  an  Inversion  by 
reciprocal  radii. 

We  apply  now  Hadamard's  formula  (6.15)  and  derive  by  comparison  of 
coefficients  near  Infinity  the  variation  of  G(P).  We  find 

c 

and  using  (12.6)  we  arrive  finally  at  tl5l 

(2.12.8)  Sc(oo)  - 4TrJJ(~gy^)2  $ ^TdC7"T 

C T 

Comparing  (12.4)  with  (12.8),  we  recognize  the  close  analogy  between  the 
functionals  C (Q)  and  C(co). 

We  can  transform  (12 .8)  by  introducing  the  tensor 

(2-12-9)  it  - «ik(7G)2  ■ a'°(p>  • 

x iC 

We  verify  easily  that 

3 ^ grji 

(2.12.9')  ).  -£**‘0  , i - 1,2,3, 

k="l  3\ 


and  we  obtain  the  formula 


(2.12.10) 


8 C(  co)  - -4TT  £ 


iii± 

n iA-1 


for  the  variation  of  the  capacity  under  a deformation  (2.2).  Thi3  result 
could  have  been  deduced  directly  from  (6.6)  and  is  applicable  for  more  general 
domains  than  is  the  case  for  (12.8) 
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Let  us  specialize  the  deformation  vector  field  and  put 


(2.12.11)  Si  “ xi  » S2  = x2 


V° 


this  deformation  leads  to  a stretching  in  the  ratio  (l+£-):l  perpendicular  to 
the  x^-axis . We  obtain  from  (12.10) 

(2.12.12)  Sc(oo)  - 8rr£  & > 

D 3 

which  proves  that  the  capacity  of  each  surface  C increases  under  such  a 
stretching.  This  result  is  not  obvious  if  the  surface  C is  not  convex,  and  it 
cannot  be  read  off  from  (12.8).  This  example  shows  the  usefulness  of  trans- 
forming variational  formulas  into  different  shapes;  for  certain  types  of 
variational  kinematics,  raonotonicity  properties  become  obvious  which  are 
otherwise  hidden. 

The  above  result  can  easily  be  extended  to  more  general  affine  trans forma- 


tions 0 Let  {(a.,  ))  be  a matrix  such  that 

3 3 


(2.12.13) 


> aik^i^k~2  4-—^  aii  ^ ^ 

i,k»l  lic  1 K i-1  11  i-1 


Then,  the  deformations  based  on  the  linear  vector  field 

s.-±. 


(2.12.13' ) 


1 w “i* 


i ™ 1>2,3, 


will  obviously  lead  to  an  increase  of  the  capacity  0. 

Let  ua  study  next  the  Neumann’s  function  for  the  exterior  D of  a closed 
surface  C.  In  the  case  of  such  an  infinite  domain,  it  is  possible  to  ask  for 
a Neumann' 3 function  with  vanishing  normal  derivatives  on  C.  By  this  require- 
ment and  the  condition  of  regularity  at  infinity,  N(P,Q)  is  uniquely  determined 
and  can  be  shown  to  be  symmetric  in  P and  Q.  It  is  tc  be  expected  that  this 
Neumann's  function  will  have  a much  simpler  variational  behavior  than  the 
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harmonic  Neumann* s functions  of  finite  domains.  Indeed,  the  formulas  (10.12) 
and  (10.15)  can  now  be  applied  immediately,  just  as  in  the  case  p^O,  We  have 


(2.12.14)  SN(Q,R)  = jjVpN(P,Q)- VpN(F,R)  8VpdC7”p 

C 

and 

(2.12.15) 


’ lllSJ  1iJiCP;Q,R)  ^ 

«0  D 1,K  x 


(P) 


dT, 


For  Q near  infinity,  we  may  develop  N(P,Q)  as  follows  In  terms  of  the 

3 


coordinates  of  Q: 


3 


1 


(2.12.16) 


nip’q)  -^rtp  • 


Here  the  C^(P)  are  harmonic  functions  of  P in  D except  at  infinity;  they 


possess  near  infinity  a development 
(2.12.17)  ^(P)  - x±  + r~3  ^ y <X,lkxk+  °(r~3) 


and  satisfy  on  the  boundary  C of  D the  conditions 

a<p,(p) 

(2.12.18)  027”“°  » P£C 


The  function  ^P^(P)  may  be  interpreted  as  the  velocity  potential  of  an 
irrotational  incompressible  fluid  flow  in  the  region  D, bounded  by  the  rigid 
wall  G,  which  has  at  infinity  unit  velocity  in  the  direction  x^.  The  coeffi- 
cient matrix  ((CX^))  can  be  shown  to  represent  a tensor;  it  is  customary  in 
hydrodynamics  to  consider  the  tensor  C 19 1 


(2.12.19) 


w.,  = 4TTC*ik“  V 


ik 


'ik 


where  V is  the  volume  enclosed  by  G.  We  call  ((w  , ))  the  tensor  of  virtual 

ik 

mass  of  the  surface  C;  it  plays  an  important  role  in  the  hydromechanics  of 
the  surface  C moving  in  an  infinite  incompressible  fluid. 
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We  derive  from  (12.14)  and  (12.16)  by  a comparison  of  coefficients  of 

"3  • 

^i  P Var^at^ona~*-  formula  for  the  velocity  potentials 

(2.12.20)  Scpi(R)  “ JJVcj)i(P)-^7pN(P,R)Svpd<r“p  , 

C 

and,  analogously,  we  obtain  from  this  formula  and  (12.17)  by  a further 
comparison  of  coefficients 

(2.12.21)  Sc^k-ZTf  ' 

c 

These  results  could  also  have  been  derived,  of  course,  from  (12.15). 

We  would  have  found,  at  first,  the  variations  of  ^ and  the  (X^  expressed 
as  integrals  over  D and  have  deduced  (12 .20)  and  (12. 2l)  by  integration  by 
parts . Since  the  variational  formulas  involving  domain  integrals  have  an 
interest  in  themselves,  we  want  to  give  at  least  the  variational  formula 
for  the  tensor  (12.19)  of  the  virtual  mass.  We  define 


) 


(2  12  22)  ,(l.k)  _ 3<?1  9<?k,  £fEl  gl?k  5.  (CTu,  ,VCB 

(2.12.22)  ^ 9*/  9*m  9,  * WVtfiVCPk- 

which  is  for  fixed  i,k  a symmetric  tensor  in  ^ and  m.  It  is  easily  verified 
that  its  divergence  vanishes. 


(2.12.23) 


c7  21lL 


(i,k) 


m-  1 


m 


Consider  now  a variation  (2.2),  but  assume  that  the  functions  S (P)  vanish  at 

-3  ^ si 

infinity  at  least  like  r j in  this  case  the  integrals  of  — * , extended 

x 

over  the  infinite  domain  D,  converge.  Under  such  a variation,  we  may  replace 
(12.21)  by 

(2.12.24)  «7=r\ll  tr»~'  dT 


_ _s_  rfT  yz  t«i,k)  —i 

^Itc  ATrJjJ^V  3x 


m 


m 
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We  find  similarly  for  the  variation  of  the  volume  V 

(2.12.25)  $v  * -e  dT, 


JjJ 


and  consequently,  by  definition  (12.19), 


(2.12.26) 


S w 


ik 


as. 


■ <£’k)  + sik  W FT  « 


fTnKL  ' m 

V 

This  formula  is  of  value  in  the  case  of  extremum  problems  involving  the 
virtual  mass  if  it  is  not  sure  a priori  that  the  boundary  surface  C of  the 
extremal  domain  D is  regular  enough  so  that  the  variational  formula  (12 .21 ) 
applies.  Since  extremum  problems  of  this  type  play  an  important  role  in  the 
theory  of  discontinuity  surfaces  in  fluid  dynamics,  the  interest  of  the 
above  formula  is  obvious . 


If  we  study  the  development  (12.5)  of  the  Green’s  function  near  infinity 
in  more  detail,  we  are  led  to  a series 

3 

(2.12.27)  G(P,E)-^*  — 1— r 2Z  H'.Wq,  * 0(  p"3) 

P 4tt/  1-1  1 1 1 

where  the  ^(P)  are  harmonic  functions  for  P€  D,  except  at  infinity.  There, 
they  have  a development 

(2.12.28)  ^(F)  » xi  + J^L  + r“3  XI  £ik  x}c+0(r“3) 

and  they  vanish  on  the  boundary  surface  C.  They  form  the  electrostatic 
counterpart  to  the  velocity  potentials  ^(P)  derived  above  and  a completely 
dual  variational  theory  can  be  given  for  them.  This  theory  was  developed  in 
the  Hadamard  notation  in  [191  and  can  be  readily  transformed  into  the  domain 
integral  form  involving  variational  tensors . 
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13.  An 

We  shall  now  apply  (12.24)  to  a modification  of  the  extremum  problem  in 

potential  theory  discussed  in  Section  2 of  Chapter  I and  shall  illustrate 

thereby  the  value  of  interior  variational  formulas  in  such  extremum  questions. 

We  consider  a closed  surface  C which  bounds  a finite  body  B and  which  is 

o J o 

homotopic  to  a torus.  We  shall  seek  another  surface  spanned  through 

so  that  the  surface  C +C.  has  an  exterior  without  irreducible  closed  curves, 

o 1 

and  such  that  the  coefficient  pf w Ofii  connec^e^  with  the  virtual  mass  of 


C +C,  is  a minimum, 
o l 


are  not  able  to  discuss  here  the  existence  of  an 


extremal  surface,  but  we  shall  concentrate  on  the  necessary  conditions  for 
such  a minimum. 

We  define  a deformation  vector  field  S^(P)  as  follows.  We  choose  a 
point  X not  on  0^  and  describe  around  X a sphere  K^(X)  of  radius  p $ we 
suppose  p so  small  that  Kp  does  not  intersect  the  surface  C^.  We  determine 
next  a function  H(P)  which  is  twice  continuously  differentiable  in  the  entire 
space,  has  the  constant  value  1 in  the  sphere  Ky^X)  and  vanishes  outside  of 
the  sphere  . Let  Q be  an  arbitrary,  but  fixed,  point  in  not  on  C^j 

we  then  introduce  the  deformation  vector  field 


(2.13.1)  S±(P)  » 8 . J H(P)  , r = PQ  , J fixed  . 

This  field  has  a singularity  at  the  point  Q,  but  it  is  easily  seen  that  the 

variational  formula  (12.24)  still  holds,  although  the  integrals  involved  ars 

now  improper.  The  variation  (2.2)  based  on  the  vector  field  (l3.l)  will 

* 

transform  the  surface  into  a new  surface  G^+C^  which  is  a competing 

surface  with  respect  to  our  extremum  problem.  Hence,  its  corresponding 
* 

functional  (ft  cannot  be  less  than  the  (ft  of  the  surface  C^+C^,  whatever  the 
sign  of  the  parameter  £ -may  be.  Hence,  we  find  by  (12.24) 
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(2.13.2) 

and  by  (l3.l) 

(2.13.3) 


If  i 


es, 


D 


dc 


u 


D 


\ S *j.  it  ^ H(p”d'c  ■ 0 

m»  1 J m 


Because  of  the  definition  of  H(P),  the  last  integral  is  only  to  be  extended 

over  the  intersection  of  the  extremum  domain  D with  ths  sphere  K„  around  the 

r 

fixed  point  2.  If  Q£K  . as  assumed,  we  may  put  (13.3)  into  the  form 
3 / nd 

(2.13.4) 


i:  I 

m*  1 v 


9 


(izr)dtT 


where 
(2.13.5) 


$</>) 


W 2 


pq  p *1 


$‘P)(Q) 


j is  a harmonic  function  of  Q in  Clearly,  also 

SSv^(t)atp-^(Q)  ’ 3-1’2’3’ 


represents  a harmonic  function  of  Q in  far  Suppose  now  that  the  sphere 

intersects  the  surface  G^j  by  our  construction,  the  function  ^(Q)  is  regular 

harmonic  inside  even  as  the  argument  point  Q moves  across  0^.  If  Q lies 

away  from  C^+C^,  it  is  clear  that  $^(Q)  is  regular  harmonic,  since  the 

divergence  of  vanishes  identically  in  D.  Thus,  we  have  proved  the  following 

theorem:  A necessary  extremum  condition  in  our  problem  is  that  the  integral 

sums  (l3.5)  be  regular  harmonic  functions  of  Q outside  of  the  given  body  B^. 

We  may  consider  this  condition  as  a set  of  three  singular  integral  equations 

for  the  variational  tensor  t^.  »e  should  expect  intuitively  that  the  surface 

0,  is  a sheet  spanned  through  the  torus  B such  that  C +C,  encloses  the  same 
X o o 1 

volume  as  did  alone . This  follows  from  the  monotonic  dependence  of  Of  upon 
the  domain.  If  we  suppose,  moreover,  that  is  a smooth  surface,  we  recognize 

easily  that  the  extremum  conditions  on  the  functions  mean  exactly  that 

2 ^ 

(VCpx)  is  continuous  across  the  sheet  C as  was  indicated  in  Chapter  I, 
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Section  2.  Thus,  (13.5)  is  an  essential  intermediate  result  formulating  the 
decisive  extremum  condition  before  the  smoothness  of  the  extremum  surface  is 
ensured. 

It  is  quite  instructive  to  study  the  plane  problem  which  corresponds  to 

the  above  question.  We  give  two  closed  curves  and  C'  in  the  plane  and 

want  to  connect  them  by  a curve  C.  such  that  the  exterior  D of  C +0  +C 

'1  o o 1 

has  a minimum  virtual  mass,  say  in  the  x^-direction.  It  is  easily  seen  that 


the  preceding  reasoning  leads  to  the  extremum  condition 
2 

log  zb"  dt„  - $,(Q) 

PQ  J 


(2.13.6) 


£51 

m-1 


<9x 


m 


where  C^(Q)  is  regular  harmonic  in  the  common  exterior  of  and  G*f  and 


(2.13.7) 


99.  99 

3?-  (V<Pir  & 


5 


^p^(P)  being  the  velocity  potential  in  the  direction  of  the  x^-axis  for  D. 

In  order  to  solve  the  set  of  integral  equations  (13-6),  we  differentiate 
this  equation  with  respect  to  and  use  the  well-known  formulas  for  the 
interchange  of  differentiation  and  integration  in  the  improper  integrals  of 


potential  theory.  We  obtain 


(2.13.8) 


-Tit 


SSI 


jm  9x^£x 


m 


(log  _ )dt- 

PQ  r 


Observe  now  that,  by  definition  (13 - 7) - we  have 


(2.13.9) 


*11  + *22 


0 


*12  " *21 


2 1 

Using  further  the  fact  that  V log  ■ 0,  we  then  find 

PQ 


(2.13.10) 


aqi 

5^1  , l$2 

3q2  0ql 


■ TT(t22-  hi’ 


* *2i ) 


“2rrt 


11 


*2  n't 


12 
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Thus  we  succeed  in  the  plane  case  to  get  rid  of  the  improper  integral  terms  by 

considering  a proper  combination  of  derivatives  of  the  harmonic  functions 

®,(Q).  The  new  formulas  make  it  evident  that  the  tensor  components  t.,  are 
j lie 

harmonic  functions  in  the  original  given  domain  D^. 

The  formulas  can  be  simplified  still  more  by  the  use  of  complex 
notation.  We  define  the  complex- valued  function 


(2.13.11)  F(z,z)  = ^(P)  + i |>2(P) 


Then,  the  two  equations  (13. 10)  can  be  united  to  yield 


(2.13.12) 


73  3 2 

_ F(z,z)  - -'fr(t11+  it^)  - -4(~ _ ) 


d z 


using  the  complex  operators  (8.4).  Since  F(z,z)  is  a harmonic  function  and 
satisfies 


U.13.13) 


02F 

07  0Z 

d? 


0 


we  conclude  that  — 2 is  independent  of  z and  hence  is  an  analytic  function  of 
_ 3 z 

z in  D . Hence,  finally,  we  obtain  from  (13.12)  the  extremum  condition  that 


(2.13.14) 


( 


^Cp  (z)  2 


3z 


) » $(z) 


is  a regular  analytic  function  of  z in  DqJ 

The  same  method  applies  to  many  extremum  problems  of  potential  theory  in 
the  plane,  and  this  fact  explains  the  easy  application  of  interior  variations 
to  such  problems . The  analogy  and  the  difference  between  the  treatment  in  the 
plane  and  in  space  becomes  obvious.  The  same  accident  which  permits  conformal 
mapping  and  analytic  functions  to  play  a role  in  the  potential  theory  of  the 
plane  is  also  responsible  for  the  simple  results  of  the  variational  technique 


in  this  case. 
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CHAPTER  III 

THEORY  OF  THE  SECOND  VARIATION 


i. 

The  method  of  interior  variations  based  on  a deformation  field  (2,2.2) 
enables  us  to  calculate  also  higher  order  variations  of  Green’s  functions  in 
a simple  way.  In  fact,  it  is  easily  seen  that  wo  can  obtain  from  the  basic 
identity  (2.4.14)  as  many  derivatives  of  g(Q,Rj  £)  with  respect  to  the 
parameter  £ as  the  coefficient  p(x^)  in  (2.1.1)  has  Holder  continuous  partial 
derivatives  with  respect  to  the  variables  x^. 

We  shall  deal  in  this  chapter  with  the  theory  of  the  second  variation  of 
the  Green’s  function  and  derive  an  elegant  expression  for  it.  The  significance 
of  this  result  for  the  treatment  of  extremum  problems  and  for  the  general  theory 
of  the  domain  dependence  of  the  Green's  function  is  obvious.  In  this  section, 
we  prepare  the  way  by  establishing  some  useful  identities  which  have  an  interest 
in  themselves . 

We  start  from  the  identity  in 


(3.1.1) 


Lp[g(Q,R;  e)]  = o 


differentiate  it  with  respect  to  £>  and  put  £ “0.  We  obtain 


(3.1.2) 


t r 2-g[Q»R»  JL 1 
L 0 € 


E-0 


] - ”L/[g(Q,Rjo) ] 


where  L is  the  differential  operator  defined  in  (2.1.1)  and  1/  ■ if  is  defined 
by  (2.4.27)  and  has,  in  view  of  (2.6.4),  the  form 

r-"7  » 0 s 3 s a S,  - 

(3.1.3)  V[u]  - 2L.  S,k  Z-,  -xr1  - -xt1  - -rf5)  -If-) 


i,k-i  3xi 


3-i  axj  0xk  a\ 


- XI  -gf-  (pS,)u 
3-1  3 3 


This  rather  complicated  expression  simplifies  considerably  if  u is  a solution 


of  (2.1.l).  In  this  case,  we  have 

3 

(3.1.4)  L'tu]  - -LtT"'  ^ s 1 

w 3xj  J 

Since  G(Q,R) * g(Q,R;o)  is  a solution  of  (2.1.1),  we  derive  from  (3.3.2) 
the  representation 


(3.1.5) 


_g>g(Q,R:  Z ) 


9C 


s (Q)  + s (r)  ] + H(Q,H) 


a r 


where  II(Q,R)  is  a symmetric  function  of  both  argument  points  which  satisfies 
in  dependence  on  each  the  partial  differential  equation  (2.1.1).  It  can  easily 
be  seen  that  H(Q,R)  is  continuously  differentiable  even  for  Q*R,  so  that  the 
first  right-hand  term  in  (1.5)  contains  the  singularity  of  at  the  point 

R. 


Since  g(Q,Rj£)  = 0 in  £ for  QeC,  R£D,  we  have  also 

(3.1.6)  £ ? , 0 f for  Q£G,  R€D  , 

o 

and  hence  ^ 

(3.1.7)  H(Q,R)  = - (Q)  , for  Q6C,  R£D 

k- 1 ^ qk  K 


Since  H is  a solution  of  (2.1.1),  we  can  express  this  function  in  terms  of 


its  boundary  values  (1.7)  by  means  of  the  Green’s  function  and  find 


(3.1.8) 


H(Q,R) 


Comparing  this  result  with  (2.6.15),  we  recognize  the  significance  r f the 
regular  function  H(Q,R).  It  represents,  up  to  the  factor  £.  , the  first 
variation  of  the  Green's  function  G(Q,R)  under  variation  of  the  domain. 
This  function  will  play  an  important  role  in  the'- variational  theory. 
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2 • 2M.  r-e 


We  consider  the  formula  (2.4.28),  which  holds  identically  in  6,  and 
differentiate  this  identity  with  respect  to  the  parameter.  Since  the  singularity 
of  the  function  is  well  known  by  formula  (l.5\  it  is  easy  to  establish  the 

validity  of  this  process „ We  have  to  assume,  however,  that  p(x^)  is  twice 

It 

Holder  continuously  differentiable  with  respect  to  all  its  variables,  sines 
these  second  order  derivatives  will  occur  in  our  formulas.  After  differen- 


tiation, we  put  £ * 0 and  obtain 


(3.2.1) 


8 £ I £-0 


-iff*-’ 


- JJJe'[g(p5q)  , ig^£Ri 


, G(PjR) Id'C, 


E"[G(P,Q),G(P,R)!dtL 


Here,  we  put  E/ » E^  a^d  define  E^  by  (2,4,29);  analogously,  we  define 

(3.a.2)  ^Iu.,1  - II  3 W e ) 8%i  ^ « 

and  put  E*  * E ", 
o 

Consider  the  expression  ^C(P,Q)a  and  -the  identity 

£=■0 

3 0S  „ (R) 

E'C^P.Qj/KFjRjfctp  - - ~r  ^R , Q ) 

D 3 * 1 J 


which  follows  from  Green’s  identity  and  the  fact  that  ^4J(P»Q)  vanishes  on  the 
boundary  C of  D.  By  virtue  of  (1,2),  we  may  write  this  identity  also  in  the 


form 


rrr  V”5”  9 s (R) 

(3.2.4)  \ U'['V(P,Q),G(PtR)TdCp  - ■ 'ViR.Q)  2 - Z~' 

°D  " J-l 


7.(p?a)L  [%(p,R)]d'r1 
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Using  again  Green's  identity  and  the  singularity  term  of  ^ given  by  (l.5), 
we  obtain  finally 

(3.2.5)  JjJ  E'[^(FsQ),G(P?R)ld1^  = -Qo[%(P,Q),  ^F,R)] 

J D 

where  Q is  the  bilinear  form  belonging  to  the  Uirichlet  integral  (2.2.101). 
o 

The  second  term  on  the  right  side  of  (2d)  can  be  treated  in  the  same 
way,  and  we  obtain  therefore  finally 


(3 .2.6) 


3 g(Q,R:S). 


£^0 


9Q  r^.g I3LSU.&1 

0 5 ££  £ = 0 


jjJ  B"[G(P,Q),G(P,R)  Jd- 


This  important  formula  permits  an  easy  determination  of  the  second  variation 
of  the  Green’s  function. 

For  the  sake  of  completeness,  we  bj.vJ  here  also  the  expressions  for  the 


g2  3 2 

terms  — E— r A.,  and  — 2 — r P . These  expressions  are  obtained  by  straight- 
ik  2 Z,2  1 


forward  computation  from  the  definitions  (2.2.8)  and  (2.2.8').  We  find 


(3.2,7) 


1 B\*ur>  1 1 


at 


J-  3(s  ,s  ) 0(s  S ) 

■ 2-  t-55^  + •FutZT  ’ 


£“0  j * 1 


V J 


I'  J 


v-  ds.  ds  ?(  s s ) d(s  sj  3(  s s ) 

+ > — i — k+o  r_ 11_2_  + — 71-3-  4 21  3..  1 

dxj  dx.  °ik  ^(x^x^)  8(x^,x^)  ^(x^yX^)  ■ ’ 

, 92/0(x„;  E.)  3(S,,S  ) 3(S  3 ) 3(S  ,S  ) 

(3.2.7O  j jp—  - p(^)Cy(^7  * TG^i  * 


0S 


T*  [— j i ^ s + 1 
1— * L 1 ? 


E- 


S.  S j 


i, j-1  ^Xi  J 2 ^xi  ^Xj  1 J 

We  observe  next  the  particular  structure  of  the  function  g(Q,R; £ ) as 
defined  by  (2.2.3).  We  have  the  series  development 
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(3.2.8) 


g(Q5Rj  C)  = G*(Q,R)  + £ 2Z  <ir  S fQ)  + |f  S fR))+...  , 

j-l  *qi  ° i 


where  G is  the  Green* 3 function  of  the  domain  B obtained  from  D by  the 
deformation  (2.2.2).  Suppose  now  that  we  change  the  vector  field  S^(j^) 

arbitrarily  in  D,  but  keep  it  unchanged  on  the  boundary  C and  at  the  two 

* 

points  Q and  R,  Under  such  a change , D would  remain  invariant,  and  hence 

G (Q,R)  would  be  unaffected.  In  view  of  (3.2„8),  we  would  also  find  g(Q,R;£) 

preserved  under  this  change  of  the  deformation.  Hence,  we  conclude  that 

g(Q,Rj£  ) depends  only  on  the  values  of  S on  G and  at  the  points  Q,R.  The 

2 ^ 

same  holds  also  for  . Thus,  we  know  a priori  that  we  can  trans- 

ae* 

form  the  volume  integrals  in  (2.6)  into  surface  integrals  extended  over  the 

boundary  C;  under  the  process  of  integration  by  parts,  the  singular  points 

Q and  R will  yield  certain  residua  which  will  lead  to  the  special  dependence 

of  the  left-hand  term  on  these  two  points.  It  is  possible,  in  fact,  to  obtain 

such  a formula  by  starting  from  (2.3).  In  view  of  the  complicated  structure 

of  the  expressions  (2.7)  and  (2.7')  this  procedure  Is,  however,  rather  involved; 

we  proceed,  therefore,  in  a different  way. 

We  assume,  at  first,  that  the  boundary  C of  D is  sufficiently  often 

differentiable  and  that  the  vector  field  is  chosen  in  such  a way  that  it 

has  at  every  point  of  the  boundary  C the  direction  of  the  outward  normal.  Then, 

* 

the  original  domain  D will  lie  in  the  deformed  domain  D and  the  Green's 
# * 

function  G (P,Q)  of  D will  be  a solution  of  (2.1.1)  everywhere  In  D,  except 
at  the  point  P*Q.  For  P£C,  we  have  the  condition 

(3.2.9)  G*(P*,Q)  - G*(P,Q)  + E.  YZ  S (P) 

i - 1 9 pi  1 


2 u£i  0pi  apk  1 


S,(P)S,  (P)  + o( G ) - o 
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!) 


Hence,  the  function 

(3.2.10)  Q(P,Q)  - G*(PSQ)-G(P,Q) 

is  a regular  solution  of  (2.1.l)  in  D with  the  boundary  values 
for  P 6 C . 

Therefore,  we  may  represent  0(P,Q)  by  means  of  the  Green’s  function 
G(F,Q)  and  the  above  boundary  values.  Since  has  the  exterior  normal 
direction,  we  find 

(3.2.11)  Q(«.r)  - + e SAJS^l  |S(p)|  -2^^  xrT 

c p 

2 JJ  dl/2  dVf  P 

We  simplify  this  result  by  use  of  the  function  H(Q,R),  defined  in  Section  1. 
Observe  that  I Sj  ~-(S-p)  and  hence  by  ( 1 . 8)  and  C 1 - V) 

(3.2.12)  Q(Q,R)  - e H(Q,R)  + £ ^ H(P,R)dCTp 

c 


{s.v)hcr?*o(t2) 


£_  _ 

2 )J  3^ 

C P 


Since  0.(Q,R)  is,  up  to  higher  order  terms  in  £,  replaceable  by  £H(Q,R),  we 
may  write  finally 

(3.2.13)  G*(Q,R)  - G(Q,R)  = £H(Q,R)-  £2Qo[H(P,Q) ,H(P,R)  ] 

- -f  (T  J&iZfl  m&Sl  (s . tfi<r  ♦ o(  e2)  . 

2 JJ  2V2  a2/P  F 

C P 

Thus,  we  have  proved  under  the  above  restricting  assumptions  the  formula 
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(3.2.14)  S2G(Q,R)  - -2£2Q  [H(P,Q),H(P,R)1-  £2  1 dff\ 

C dV? 


2 " ~P 


32G 


In  the  last  integral,  we  may  replace  t.  by  lower  derivatives  if  we  use 

the  fact  that  G(P,Q)  satisfies  the  differential  equation  (2.1.1)  and  vanishes 
identically  on  C.  If  p^  and  denote  the  principal  radii  of  curvature  at 
the  point  P,  it  Is  easily  checked  that  in  view  of  (2.1.1) 


(3.2.15) 


G.M ZAL  „ f_i_  + JL)  &MZA1 

\ n /o  / A * 


A Pi 


Hence,  we  arrive  finally  at  the  result 

(3.2.16) 


82G(Q,R)  - -2£2Q  [H(P,Q),H(P,R)1 
o 


P2  f f gotP.tt)  JlaiZiEl  i- 1_  ♦ -L.us.i/)2i/r 

8 jJ  3A  9Vf  '/^  />2,ls  d<rP  • 


We  want  next  to  rid  ourselves  of  the  assumption  that  the  deformation 

vector  field  has  everywhere  normal  direction  with  respect  to  the  boundary  C. 

Let  us  suppose  that  only  the  assumption  (-^*S)-<0  is  fulfilled  on  C.  It  is 

easy  to  determine  the  normal  shift  A on  C which  is  necessary  in  order  to 

transform  this  surface  into  the  surface  C induced  by  the  deformation  (2.2.2). 

In  fact,  let  x^(u,v)  be  a parametric  representation  of  the  surface  C and  let 

Z^(u,v)  be  the  components  of  the  unit  vector  on  C in  the  direction  of  the 

# * 

interior  normal.  We  want  to  represent  the  point  x^  of  C in  terras  of  a 
normally  displaced  point  x^(u+  6 V)  on  C.  Thus,  we  have  for  a suitable 

X*  * 

of  u « u+  £ U and  v =>  v + £ V the  condition 

$ X 

(3.2.17)  xi(u,v)  + £ S^tujv)  =>  x^Ujv)  + £ (”^~  U + V) 


,.2  ^2x.  , ^2x 

+ r --4  ]tt  + ? — 

2 L ^ 2 0u  £v 


L. 


~2 
d x 


UV  + 


- V21 


'd'V . 


eA(u*,v*)  -^(u,*)  + e u + v )/■■*■  °(  &2)  * 
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Multiply  this  equation  with  ~V,  (u,v)  and  sum  over  i.  Since  U is  orthogonal 

d 9 x^  3 2'^  9 ^ 1 

to  the  tangential  vectors  — , “^T,  ~9u  5 *~9v  ’ we  obtain 

(3.2.18)  A(u*,v*)  = (/U2^  2?HUV+7tV2)  + o(£  ) , 

where  ^ 7\  are  the  coefficients  of  the  second  fundamental  form  of  C.  We 

deduce  from  (2.17)  by  comparison  of  the  tangential  components 

(3.2.19)  UE  + VF  - (S*  -^2)  + o(l)  , UF+ VG  » (S-  ~^)  + o(l) 

o to  £7  V 

where  E,F,G  are  the  components  of  the  first  fundamental  form.  Thus,  U and  V 
can  be  expressed  linearly  in  terras  of  S and  geometric  quantities  of  C,  up  to 
higher  order  terms  in  E . We  obtain  finally 

(3.2.20)  A(u,v)  - (s»  v)~  e [u  + Y 

o u O v 

- Y [/ u2  + 2WJV  71V2]  + °( e ) 

We  write  for  short 


(3.2.20')  - (5-^)-e|[s]  + o(e) 

where  (£[S]  is  a simple  expression  of  second  order  in  and 
on  the  geometry  of  the  surface  C. 

Consider  now  the  expression 

(3  .2.21)  G*(Q,R)  * G(Q,R)  + 8 G(Q,R)  + j 52G(Q,R)  + o(  £ 2) 


as 

~ — ^ which  depends 
3xk 


We  have  in  formulas  (2.6.15)  and  (2.16)  simple  expressions  for  the  first  and 
second  variation  if  a normal  shift  =*  £ (S'|/)  is  made.  Inserting  the  value 
(2.20')  into  these  formulas  and  collecting  the  terms  of  order  & , we  obtain 
finally 

(3.2.22)  82G(Q,R)  » -2  &2Q  [H(P,Q)  ^(P^R)  1 


2 17-1-4. 

JJ  a^?  bpv  LVi 

G 


)(s*z/}2~  2§[s]]d<r 
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This  variational  result  has  been  derived  under  the  assumption  that 

and  that  the  surface  C is  sufficiently  smooth.  We  know,  however, 
from  our  general  considerations  that  a formula  of  this  type  can  be  obtained  by 
integration  by  parts  of  the  expression  (2.6).  We  used  the  later  reasoning  only 
in  order  to  save  laborious  calculations  and  to  arrive  easily  at  the  final  result 
of  these  transformations.  Hence,  formula  (2.22)  must  hold  for  an  arbitrary 
deformation  vector  field  and  for  all  surfaces  C which  possess  continuous 
curvature  and  a continuous  second  normal  derivative  of  the  Green's  function. 

In  particular,  the  elegant  formula 


(3.2.16’)  S2G(Q,R)  - -2Q  [S  G(P,Q) , £g(P,R)  1 

o 


^G(P.Q)  ^G(P,R)  (-L-  + 
dv  ^ p1 


dcr 


holds  for  all  such  surfaces  under  a normal  shift  . 


We  calculated  in  the  preceding  section  the  second  variation  of  the  Green's 
function  of  a domain  D under  a deformation  of  the  particular  type  (2.2.2).  It 
Is  necessary  to  adapt  this  variational  formula  to  various  problems  arising  in 
the  applications.  For  example,  we  may  have  a family  of  closed  surfaces 
which  depend  on  a parameter  t;  they  enclose  domains  with  corresponding  Green's 
functions  G (Q,R).  It  is  necessary  to  find  simple  expressions  for  the  first  and 

Xr 

second  derivatives  of  G^_  with  respect  to  the  parameter  t.  These  expressions  are 
obtained  from  our  variational  formulas  by  elementary  transformations . 

Let  us  suppose  that  two  surfaces  Gq  and  are  given  such  that  encloses 
Cq  entirely,  and  let  V(P)  be  a twice  continuously  dlfferentiiible  function  3n 
the  domain  bounded  by  and  Gq0  We  assume,  moreover,  that  |VV|/o  in  this 
region  and  that  the  function  V(P)  has  on  and  the  boundary  values  0 and  1, 
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respectively.  Consider  now  the  level  surface  C , defined  as  the  locus 


(3.3.1) 


V(F)  = t , Pec, 


These  surfaces  have  no  self-intersection  and  they  determine  a one-parameter 

family  of  finite  domains  which  depend  monotonically  upon  the  parameter  t. 

Let  G^(Q,R)  be  the  corresponding  Green's  function;  our  problem  is  to  determine 

the  derivatives  of  G with  respect  to  the  parameter  t. 

t 

Let  P be  a point  of  C and  compute  the  normal  shift  $2/  which  is  necessary 
in  order  to  place  P on  the  level  surface  We  clearly  have  the  condition 


(3.3.2) 

whence 

(3.3.2') 


S*'  + 2 ~ — 2 ( § ^ )2  + . • . * A* 


S*.  (-Hr)-1  At  - i At2*°CAt2) 


Therefore,  we  derive  from  (2.6.15)  and  (2. 16')  the  identities 

a CC  ^C.(PjQ)  3G  (PSR)  ov/v)  i 

(3.3.3)  -2-  a (n  v)  - . U t i— 


-fi  Gt(4’R! 


v 9*  ’ dcrp 


(3.3.4) 


fs  Gt(,j>R)  • Gt(p’Q)  > -ft  Gt(p-R);i 


0G+( P,Q)  aG+(P,R) 


_Lj 

3V  A /°2 


" • <f£rl»rP  . 

We  observe  now  that  C is  a level  surface  for  the  function  V(P).  Hence,  we 

t 

have  the  well-known  identity 


(3.3.5) 


v2v  = -^-  Hr 


* _av_ 

a/ 


and  we  may  write 
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(3.3.4') 


2 G+(Q,R)  - -2Q  1-%Z  G (P,Q),  ^7  G (P,R)] 

l *D  O y u u ^ o u 


9g  (P,Q)  0G  (P,R)  2„(v)  , 2 

i 3 V V(5  >dcr 

v 9^  ; Ui.aop  . 


The  preceding  formulas  are  partlouiury  interesting  in  the  case  of  a posi- 
tive coefficient  p(x  ) in  the  differential  equation  (2.1.l).  In  this  eass,  we 

1 9Gt 

know  that  the  Green's  function  is  positive  in  the  domain  and  that  is 

9 V 

positive  on  Moreover,  it  is  clear  that  "gpT  is  negative  on  each  level 

surface  C Hence,  we  see  that  the  Green's  functions  G (Q,R)  increase 
raonotonically  with  t by  virtue  of  (3.3),  a fact  which  follows  also  easily 
from  the  minimum  principle  valid  for  the  Green’s  function.  Let  us  suppose 
further  that  the  ^unction  V(P)  Is  subharmonic,  i.e„, 


(3.3.6) 


V V 5,0 


between  and  C^.  7?e  can  then  assert  that  the  surface  integral  in  (3.41) 

is  negative  and  we  have  the  estimate 

(3.3.7)  '2L~S  G . (Q,R)^-2Q  1-2-  G (P,R)  ,JL  G (P,R) ] 

- +.  0 t 0+  + 


9t  t 


Equality  will  hold  in  (3.7)  in  the  case  of  a harmonic  function  V(P)„ 

In  order  to  draw  conclusions  from  (3.7),  we  introduce  an  arbitrary 
fundamental  singularity  S(Q,R)  for  the  differential  equation  (2.1.l)  which 
is  defined  in  the  larger  domain  bounded  by  C^.  Thus,  the  functions 


(3.3.8) 


ht(Q,R)  - Gt(Q,R)-  S(Q,R) 


will  be  regular  solutions  of  (2.1.1)  throughout  D . We  might  choose,  in 

particular,  S=*G^.  Clearly,  we  have 

2 2 _ 

(3.3,9)  ^5  ht(Q,E)  - 2-J  Ot(Q,HU  -2Q0C^  0t(P.9)  . fi  • 
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Gonsider  no?)  the  quadratic  form 

N 

(3.3.10)  XT  WV  Xf>  A<r  " kt 

p, cr-1  " 1 

based  on  N arbitrary  numbers  Xp  and  N points  Qp  inside  We  conclude 

from  (3  .9)  and  the  fact  that  Q Eu]  ? 0 for  p(x. ) ^0 

o l 

a2 

(3.3.11)  ““  A-£0  , 

9td  * 

which  shows  the  concave  dependence  of  the  quadratic  form  A^,  upon  the 
parameter  t. 

This  result  gives  an  important  and  useful  insight  into  the  dependence 
of  the  Green’s  function  upon  the  domain  D.  Consider,  for  example,  the  case 
of  Laplace's  equation.  Here,  a fundamental  singularity  is  readily  available, 
namely 

(3.3.12)  s<e>H>  * ^ irfer 

We  have 

(3.3.23)  ht(Q,Q)  « - ^ C^CQ) 

where  C (Q)  is  the  capacity  constant  of  the  domain  D with  respect  to  the 
X/ 

point  Q,  as  defined  in  Section  12  of  Chapter  II.  We  may  therefore  state  that 
in  the  case  of  a subharmonic  function  V(P)  the  capacity  constant  C^(Q)  is  a 
monotonic  and  convex  function  of  the  parameter  t. 

Let  us  return  to  the  case  p(x^)  0.  In  this  case  it  is  well  known  (ls  2] 

that  there  exists  a kernel  function  K (P,Q)  for  the  domain  D^,  and  the 
differential  equation  (2.1.l)  which  has  the  reproducing  property 

(3.3.14)  u(Q)  - Q [K  (P,Q),u(P)l 

o t» 

for  every  regular  solution  U(P)  of  (2,l.l)  in  D^.  The  Dirichlet  integral 

Q is,  of  course,  to  be  taken  over  the  domain  D . We  have  by  the  Schwarz 
o 7 x 
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inequality  and  the  reproducing  property  of  the  kernel  the  inequality 


(3.3.15) 


u(QnSQ  iu]*K,  (q,q) 
o t. 

_2_ 


Observe  now  that  ~r  G (P,Q)  is  a regular  solution  of  (2.1.1)  in  D , so  that 

<2  b t T> 

the  above  inequality  can  be  applied  to  it.  Thus,  we  deduce  from  (3.9)  the 
result  „ 

9h.(Q,Q)  2 , 

*-r-  * -2 ( Srr ) CKjQ.Q)]  1 


(3.3.16) 


9t 


We  may  define  [K  (Q,Q)]  ^ by  virtue  of  (3.15)  as  the  minimum  value  of 

Tj 

the  Dirichlet  integral  Q [u]  for  all  solutions  u(P)  of  (2.1.1)  in  L.  which 

o t 

have  at  Q the  value  1.  Hence,  this  functional  increases  monotonically  with 
Thus,  (3.16)  implies 


(3.3.17) 


a\(Q,Q)  9h.(Q,Q)  9 , 

*-5—  ^ -2( ) [K  (Q,Q)]_;L 

at  * 0 


an  inequality  involving  only  the  kernel  function  of  the  initial  domain  Dq. 
Analogous  estimates  can  be  obtained,  of  course,  for  the  more  general 
expressions  (3.10). 

It  is  possible  to  extend  the  preceding  reasoning  to  the  case  of  Laplace's 
equation.  In  particular,  if  we  consider  the  exterior  D of  a closed  surface  C, 
we  can  introduce  the  kernel  function 


(3.3-18)  K(P,Q)  - N(P,Q)-  G(PSQ)  - n(P,Q) + g(P,Q) 

where 

(3.3.19)  H(P.«)  • * »(P.Q)  , 0(P,Q)  ■ - g(P.«) 

are  the  Neumann’s  and  the  Green '3  functions  of  D. 

Let  Cq  be  a closed  surface  including  0^  and  let  V(P)  be  subharmonic  in  the 
domain  between  the  two  surfaces.  We  assume  again  that  iWl/O,  and  that 
V(P)»0  on  C and  V(P)  ■ l on  C_ . Let  D be  the  exterior  of  the  level  surface 
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V=»t  and  let  G^_  be  the  corresponding  Green’s  function.  We  find  in  analogy  to 
(3.9)  the  inequality 


(3.3.20) 


320t(oo) 

0t2 


^ 8TTQ  [ 
o 


-2. 

3t 


Gt(P)] 


for  the  capacity  G of  the  domain  D , where  G (P)  is  the  Green's  function  of 

b U 

D.  with  the  singularity  at  infinity.  Since  ~~  G.  (F)  i3  regular  harmonic 
u yt  t 

in  D^,  we  have 

with  K the  kernel  function  of  D and,  hence,  by  the  Schwarz  inequality 

(3.3.22)  [*df;  G (Q)]2-£K  (Q,Q)Q  {-fr  G (P)] 

(7t  t t O d Z t 

Combining  (3.20)  and  (3.22),  we  find 

02C  ( oo ) p 

(3.3.23)  ^8Tr[^t  Gt(Q)fKt[Q,Q]"i 

<2 1 


for  arbitrary 
(3.3.24) 


Q e . Let  now  Q 

lira  r2K  (Q,Q)  =» 
Q— >co 


co;  it  can  be  shown  that  [19  ,p.  139 1 


Ct(  00) 

4TT 


while  obviously 

(3-3.25)  lim  r -fj  0 (Q)  • - ^ 

y ~>cd 


3Ct(co) 

at 


Hence,  we  find 
(3.3.26) 


a2ct(p) 

at2 


^2Ct(oo) 


, 3C  ( 00)  p 

•If 1 ^2 


at 


or 

(3.3.27) 


9 


Thus  the  reciprocal  of  the  capacity  varies  as  a concave  function  under 
a level  surface  variation  with  a subharmonic  function  V(P)„  This  result  shows 
the  significance  and  usefulness  of  the  variational  formula  (3.4*). 
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We  can  use  (3-27)  in  discussing  various  extremum  problems.  Let,  for 
example,  C be  a closed  surface  which  encloses  a fixed  point  Q.  Define  the 
functional 

(3.3.28)  TT(C)  - C(Q)-C(co)  , 

which  is  dimens ionles 3 , i.e.?  remains  unchanged  under  a similarity  transfor- 
mation of  the  surface  C. 

We  call  a smooth  surface  C stationary  with  respect  to  the  functional  YT 
if  we  have  ■ 0 under  infinitesimal  deformations  Sz'of  C.  By  the  Hadamard 
formulas  (2.12.4)  and  (2.12.8)  this  condition  has  the  form 

(3.3.29)  C(Q)  2 — 2-  G(P,Q)  - C(oo)  2 G(P)  , 

dP~  3p+ 

» + 

for  every  Pe  C,  where  P and  p denote  the  interior  and  exterior  normals 
of  C,  respectively. 

If  C is  a sphere  of  radius  R around  Q,  we  have 

(3.3.30)  G(P,Q)  - “ £)  , G(P)  - (1-  J 

hence  ^ 

(3.3.31)  C(Q)  2 G(P,Q)  » Kj5^C(cd)  2 ^rG(P)  . 

4TTR- ' VP 

Thus,  the  sphere  is  a stationary  surface  with  respect  to  the  functional  IT . 

We  want  now  to  show  that  it  is  the  only  possible  one . 

Suppose,  indeed,  there  were  another  surface  with  this  property.  We 

take  to  be  a sphere  inside  with  sufficiently  small  radius  and  we 

consider  the  level  surfaces  0^(0^- 1 ^a.  l)  of  the  harmonic  function  V(P)  which 

is  zero  on  C and  1 on  C. . Let 
o 1 

(3.3.32)  TTt  - TT(Ct)  - ct(Q)*ct(oo) 

We  have 

(3.3.33)  TT't  = Ct(co)C'(Q)  + C+'(co)Ct(q) 
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(3  .3  = 3 4)  TTt"-  Ct(co)G^(Q)-2C^(co)Ct,(Q)  + G'/(a))Gt(Q5 

Whenever  TT^  vanishes,  we  can  write 

Cf(Q)  cr(oo)  0/(00) 

(3.3.35)  TT  - TtfjJgy  * ^ - 2<fe>2l  • 

t t t 

In  view  of  the  convex  dependence  of  C+(Q)  on  t and  by  (3-26),  we  find 

(3.3.36)  TT"^0  > if  TTt'  - ° 


We  have  to  consider  the  case  of  equality  in  (3.36).  Clearly,  this  is 
only  possible  if  C^(Q)  ■ 0,  and  by  virtue  of  (3-7)  this  is  only  possible  if 

(3-3.37)  G^(P,Q)  * const.,  for  P inside  C^_  . 

By  (3.3),  this  is  only  the  case  if  we  have  on  C the  equality 

(3.3.38)  -jfr  ot(p,Q)  - , psct  , 

with  a suitable  constant  factor  . The  surface  C is  an  analytic  surface, 

V 

and  by  the  Cauchy- Kowa lew ski  uniqueness  theorem  for  solutions  of  partial 
differential  equations  with  given  initial  data,  we  obtain 
(3.3-39)  G^(F,Q)  * c*[v(p)-t] 

The  sphere  C appears  as  a level  surface  of  G (P,Q)  and  hence,  obviously,  all 
level  surfaces  C^_  are  spheres  around  the  point  Q„  Hence  the  final  surface  C^ 
will  be  a sphere  around  Q and  we  have  the  original  type  of  stationary  surface 
with  respect  to  TT. 

If  C is  not  a sphere  with  Q as  center,  we  have  the  inequality 

(3.3.40)  tt"  p-o  , if  rr'-o  . 

This  result  shows  that  IT*  can  only  cross  from  negative  values  to  positive 
values  as  t increases . Since  TT  * 0,  we  recognize  that  TT.  remains  positive 

O 0 

for  t>0,  and  TT^»0  is  impossible.  Thus,  there  cannot  be  a stationary 
surface  other  than  a sphere  around  Q„ 
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Since  TT  in-'raai*5*^  wh^n  we  passed  from  C to  G , we  see  that  IT  has  its 
t o J. 

minimum  value  for  this  spt.°'*e  Qila 
(3.3.41)  C(o,)G(Q)^1  . 


with  equality  bolding  only  in  the  cab  of  a sphere  around  Qe 

We  observe  how  the  convexity  of  the  capacity  leads  to  the  inequality 

(3.3.4l)j  a uniqueness  theorem  for  the  function"1.  oquation>A3 .29)  and  the 

\ 

monotonicity  of  TTj.  at  the  satv  time . 


4 . IheJiiq- d nal^qaas.  • 

All  methods  applied  in  the  preceding  '-actions  can  be  carried  over  wiV°ut 
any  change  to  the  case  of  a partial  differential  equation  (2.8.1)  in  two 
independent  variabl.es..  We  devote  to  this  case  a special  section  only  for  the 
reason  that  some  formulas  simplify  considerably  and  lend  uheraselves  to 
interesting  applications . 

Let  again 


(3.4.1) 


H(Q,R)  - 


aoinfil  .MlMl  (s.*,)d„ 


2m, 


dV, 


denote  the  solution  of  the  original  differential  equation  which  is,  up  to  the 
factor  £,  the  first  variation  of  the  Green’s  function  for  the  displacement 
vector  field  3^  (i-1,2).  If  S,  has  the  normal  direction  on  the  boundary 
curve  C of  the  plane  domain  D and  if  K(s)  denotes  the  curvature  of  C at  the 
point  P(s),  we  have  in  analogy  to  (2.16) 


(3.4.2) 


82G(Q,R)  - -2e2Q  [H(P,Q),H(P,R)] 


..  e2  f k(3)(s^)2 


?>v. 


ds 


We  want  to  extend  this  formula  to  the  case  when  3 is  not  normal  to  0; 
it  is  at  this  stage  that  the  two* -dimensional  theory  becomes  simpler,  in  view 
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of  the  simpler  differential  geometry.  We  decompose  the  displacement  vector 
into  normal  and  tangential  components  and  T^, 

(3.4.3)  S.  - N.+T  - (S*a/)3^  +[S  - 3 

i x i i i i 

and  we  let  N and  T be  the  projections  of  the  vector  on  the  normal  and 
tangential  vectors,  respectively.  Starting  again  from  the  condition 

(3.4.4)  x^s)  + £ s^tx^)  a x (a  )+  £ A(s  ) 24  (s  ) , 

we  derive  by  frenet’s  formulas  the  relation 

(3.4.5)  £(N,  + T^)  = ^(3)^(3  - a)-  £A(s  )(s  - s)J" 

4 '74(s)’f-^  (s  -s)2  + £ A(s  ij*  + o(  E,2) 

Thus,  we  find 

(3.4.6)  s-sss£T+o(£‘)  , 

0.4.7)  A(s*)  -N-e  “T2+o(£2) 

and  filially 

(3.4.8)  A(s)  - N-  E [T  4 -■  T21+  o(£2) 

a3  Z 

This  is  the  equivalent  normal  displacement  cf  C under  the  general  deforma" 
tion  (4.3). 

Because  of  the  correction  in  the  normal  shift,  the  first  order  variation 
leads  to  a correction  term  in  the  second  variation  and  we  find 

0.4.9)  S2G(Q,R)  - “2  62Q  IK(P,Q),H(P,R)] 

o 

" J ds 

C * r 

This  formula  gives  the  second  variation  of  the  Green’s  function  under  an 
arbitrary  displacement  with  a vector  fxeld  S .. 
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Formula  (4.9)  enables  us  to  study  the  variation  of  the  Green's  function 
Tilth  the  parameter  of  a one- parameter  family  of  curves.  Let,  for  example,  C 
be  a family  of  closed  curves  given  by  the  complex  parametric  representation 


(3.4.10)  z * f(<7";t)  , OsSt^l,  0<CT^1  , 

which  expresses  for  fixed  t the  points  of  the  curve  C in  dependence  on  the 
parameter  C T.  We  assume  the  function  f to  be  periodic  in  ( T with  period  1 
and  to  have  two  continuous  derivatives  in  each  variable. 

Each  curve  defines  a finite  domain  D^_  and  a corresponding  Green's 
function  G(z,X;t).  We  consider  the  curve  C and  the  shift 

Ti 

0 

(3.4.11)  f(<T;t)-  f(CT;t  ) - f.  (CTjt  )(t-  t ) + J f,  . (<Tjt  )(t- t )2+  ...  , 

o t o o 2 tt  7 o o 

which  carries  points  from  C into  points  of  C . We  may  express  the  unit 

u u 

O 

tangent  vector  at  the  point  z(CT*t ) in  the  form 

(3.4.12)  eiW=  f_(<T  jt  )/\Zr{<T  jt  )| 

Q o u o 

and  decompose  ? into  tangential  and  normal  components, 

(3.4.12')  f.(cr  Jt  )(t- t ) - TeiW+  Hi  ei<?<  , 

too 


with 

(3.4.12") 


O' 


Re 


ft&T  ito)e“i0t } (t-  to)  , N - Im|ft(cr;to)e"i<?fJ(t-  tQ) 


Using  (3-4*9)  and  observing  the  additional  correction  from  the  first  variation 

due  to  the  term  “ f, . (< T jt  )(t  - t ) in  the  displacement,  we  find 

2 tt  o 0 


(3.4.13) 


~ a -2Q  [H(w,z;t),H(w,^r  ;t) 


1 «i  v* ;t)|2- 2 


f Ira -ff_Jcr  ;t)e  ]ds, 


tt 


J 


'c r 
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This  formula  becomes  particularly  interesting  in  the  case  of  convex  curves , 
for  which  the  curvature  K 1b  always  positive.  Let  C and  C be  any  two  convex 

C . X 

curves  enclosing  the  origin.  We  suppose  them  described  analytically  by  the 
values  po(  <?),  Px(?)  of  the  supporting  functions  in  dependence  on  the  angle 
of  inclination  cp  of  the  normal.  We  introduce  then  the  family  of  convex  curves 
C by  defining  their  supporting  functions  as  follows: 

(3.4.14)  pt(CjO  - (l~  t)po(<p  J + tp^Cp)  , 0«£t-£.l 

In  order  to  describe  this  family  in  the  complex  form  (4-10),  we  observe  that  a 
tangent  to  whose  normal  has  the  angle  of  inclination  cp  is  given  by  Hesse's 
normal  form 

(3.4.15)  x cos  cp  + y sinCp  » p^Qp) 

and  that  C is  the  envelope  of  this  family  of  straight  lines.  Hence,  the 
points  of  C are  found  by  eliminating  Cp  from  (4.15)  and 

(3 .4.16)  -x  sin  Cp  + y cos  Cp  » p'(cp  ) 

Thus,  v;e  find  easily 

(3.4.17)  f(cpjt)  - e^Cp  (<p)  + ip'(cp)] 

- (1-  t)f  (Cp  ;0)  + tf(cp  jl) 

We  calculate 

(3.4.18)  f = iei£?Ept(cp)  + p"(cp)] 
and  verify  that  e1^  -ie1^5  and 

(3.4.19)  Im^te 
Moreover,  we  have  clearly 


(3. 4. 2.0) 


y?i‘)so 
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Thus,  (4*13)  simplifies  to 

2 

(3.4.21)  - -2Q 

- T - Gfa«a»k?  K|f(cp,-i)  - r(cp  ;0)|2  ds 

J & 2/^  C*  Z^w  Cp 

c 

Formula  (4-21)  is  very  useful  in  the  case  that  the  coefficient  p(x^)  of  the 
differential  equation  is  non-negative.  In  this  case,  we  can  again  assert  that 

O n 

Qq[u]  is  non-negative  and  that  is  non-negative.  Let  S(z,<^)  be  a funda- 

mental singularity  of  the  differential  equation  considered;  then 


(3.4.22) 


h(z,£  ;t)  =*  G( z,X  $t)-S(z,^  ) 


is  a regular  solution  of  the  differential  equation  in  D^.  Reasoning  as  in 


Section  3,  we  can  assert  that  the  expression 

N 


(3.4.23) 


At  * rtfii  h(V’ V°  Xf>  V 


is  a concave  function  of  the  parameter  t. 

This  general  result  contains  numerous  special  cases  and  leads  to  various 
inequalities.  It  implies,  in  particular,  that  the  capacity  of  a convex  family 
(4.I4)  at  a given  point  v?ith  respect  to  Laplace's  equation  is  a convex  function 
of  the  parameter  t of  the  family.  We  do  not  enter  here  into  a discussion  of 
the  various  results  and  their  possible  generalizations . 


5 . Singular  variations  of  the  second  order . 


In  Section  9 of  Chapter  II,  we  derived  from  the  general  formula  for  first 


order  variations  of  the  Green's  function  of  a domain  D in  the  complex  z-plane 


a particular  result  for  the  deformation 

* i<*  2 

(3.5.l)  z » z + ~ — —ft— 

z - z 

o 

This  formula  has  in  the  case  of  Laplace’s  equation  the  simple  form 
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(3,5.2)  G*(f  *,  fp-  G(  f 5 - Re^BTTyo 


2 ii* 
e 


<9G(z^,  £)  9u(zo,)!7  ) 


0 a 


2* 


+ o(p2) 

Since  this  result  has  been  derived  from  the  formula  (2.8.7),  it  is  valid  for 
the  most  general  plane  domains  and  can  be  applied  to  extremum  problems  for  the 
Green's  function  in  order  to  characterize  the  txtremal  domain.  One  obtains 
In  many  cases,  from  the  extremum  requirement  and  the  variational  formula, 
necessary  conditions  for  the  extremum  domain  in  the  form  of  a differential 
equation  for  its  Green's  function.  One  can  then  show  easily  that  the  required 
extremum  domain  has  a piecewise  analytic  boundary  curve  Cj  one  can  even  find  in 
many  cases  domains  D which  satisfy  all  necessary  conditions  imposed  by  the 
study  of  the  first  order  variation.  But  the  question  arises  to  show  that  there 
exists  only  one  unique  domain  satisfying  the  necessary  conditions.  It  is  natural 
to  study  the  theory  of  the  second  variation  in  order  to  obtain  such  uniqueness 
results . 

We  want  to  derive  in  this  section  a theory  of  the  second  variation  of  the 

harmonic  Green's  function  under  interior  variations  of  the  type  (5.1).  Since 

we  cannot  use  the  principle  of  superposition  of  variations  when  studying 

second  order  terms,  we  shall  start  immediately  with  a variation 

N 

(3.5.3)  z*  = z + £ H , UJ^l 

j z~ 


We  suppose  that  no  lies  on  the  boundary  G of  the  original  domain  D;  if  E 
is  small  enough,  the  mapping  (5.3)  will  carry  G into  a curve  C in  a one- to- 
one  manner,  and  let  D be  the  interior  of  C . We  denote  by  G ( ^ , tj  ) the 
Green’s  function  of  the  varied  domain. 

We  assume  that  C i3  an  analytic  curve;  by  our  preceding  remarks  it  is 
clear  that  this  case  will  be  the  most  Important  one  in  the  applications.  More- 
over, because  of  the  continuous  dependence  of  the  Green's  function  on  the 
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domain,  our  final  result  can  be  extended  to  the  most  general  domains  through 
approximation  by  analytically  bounded  ones.  We  shall  follow  the  same  method 
used  in  Section  2 for  the  general  case,  but  shall  utilize  fully  the  advantages 
of  complex  notation. 

Since  D and  C have  analytic  boundaries,  we  can  choose  £ so  small  that 

* * 
both  Green’ 3 functions  G and  G are  analytic  in  the  closure  of  D+D  . 

Writing  z =•  z + £ Z(z),  we  have  by  Taylor's  theorem 

(3.5.4)  0 « G*(£  + ez(f ),  q)  - G*(£,rj)  + 2 Rejs-Z(£ 

z(£)2  •22£--^-3-)7  + o(£3) 

J 

for  tj£C.  Hence,  the  function 

(3.5.5)  J\_(£,  7)  - G*(£,  0)-G(f,q) 


d'S 


is  regular  harmonic  in  the  closure  of  U and  has  for  ^ £ C the  boundary  values 

(3.5.4’)  -2  Re-^6  Z(f)  Z(£)2  + Q(£3)  . 

Using  the  Green’s  function  G( £ , n ) of  D and  the  identity 


(3.5.6) 
we  obtain 


JMzjJL).  d„  = i _2fi.(a,4£l 
anz  a i ^ 


dz 


(3.5.7)  A(  £,  n)  - Re-ff  \ [£  Z(z)  3,Q.ls,rii  + -|-  z(z)2  )-] 

! ! ' ^ dz 


* dz  1 + 0(  £3) 

I 


Put  now 


(3.5.7:)  A(^,  7)  - ? A^,  7)4‘J  A^f  ,q)^0(£3) 


We  find  at  once 
(3.5.8)  A x ( £ , fj) 


<9  2.  9 2 


? 
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C3Z  az 


evaluate  the  integrals  by  means  of  the  residue  theorem.  We  find, 


first,  N 

(3.5.10)  AA$,n)  - Re4  8TT  Ca, 

11  V*  1 


3G(zv,  <7)  aS(Zj,,£) 


- 2ZI 


This  result  is  equivalent  to  (5-2)  and  could  have  been  obtained  from  this 
formula  by  the  principle  of  superposition  of  first  order  variations. 

In  order  to  calculate  the  second  variation  term  )>  it  ia  useful 

to  introduce  the  two  kernels 


(3.5.11) 


(3.5.12) 


K(£,^) 


L(r  n)  . , ? G LZ*  Ql 1 M,n) 

L(r-r  4 <?*<?-?  W%£- rj  )2  * " 


which  play  an  important  role  in  the  theory  of  conformal  mapping  and  orthogonal 
function  systems  in  the  domain  D [ 1 1.  ^he  kernels  K(  ^ , rj  ) and  '"j)  ar© 

regular  analytic  functions  of  their  arguments  throughout  the  domain  D. 

In  order  to  simplify  the  result,  let  us  suppose  that  the  parameters 
in  (5.3)  have  been  chosen  in  such  a way  that 

(3.5.13)  Z($  ) - Z(r]  ) - 0 , 

i.e.,  so  that  the  argument  points  considered  are  fixed  under  the  variation. 
With  these  notations  and  assumptions,  we  calculate  easily 

r 2 9 \ \ 7}  jj j,  a,  \ 

(3.5.14)  AM,  n)  * Re 4 16 tt  l___^  Avy  — — 2 

2 Y 1 L h.^-1  v i* 


16TT2  S , A v A 


9G(z u,n)  &G(zu,  $ ) 


M "vy  'biy 

yV™  l * 


K(*y  zp) 


JN 

* Sir  Ay  A 


2 G (2^,77  ) 0G(zy,  >^)^  3u(zp,  ^ ) ^G(z^> 


* ' 1 V M / y 

\xy  **I  Zy) 
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In  this  formula  for  the  second  interior  variation  of  the  Green's  function, 

certain  bilinear  forms  appear  which  are  well  known  in  the  general  theory  of 

conformal  mapping.  The  forms 

N 

(3.5.15)  L«  / ^ i(z  tzp)x  , K - 

u,-p  -l  r ' 

satisfy,  for  example,  the  inequality 

(3.5-15’)  Il'I  ^.K 

a result  which  leads  to  various  convexity  statements  with  respect  to  functionals 
connected  with  the  Green's  function. 

We  do  not  intend  to  make  applications  of  the  above  formulas  in  the 
present  paper.  Our  main  purpose  in  presenting  the  results  is  to  show  the 
uniform  character  of  the  various  methods  applied  in  the  variational  calculus. 

In  this  section,  we  want  to  give  an  alternate  method  for  the  derivation  of 
second  order  variational  formulas  and  apply  this  method  to  the  Neumann's 
function. 

We  deal  again  with  a three-dimensional  domain  D and  the  partial  differential 
, . * 

equation  (2.1.1).  Let  D be  another  domain  and  suppose  that  the  intersection 

* D 0 D is  a domain.  Let  N(P,Q)  and  N (F,Q)  denote  the  Neumann’s  functions 

■59 

of  the  domains  D and  I)  , respectively . We  continue  the  definition  of  these 
functions  beyend  their  domain  by  putting  them  equal  to  zero  if  an  argument  point 
lies  outside  of  this  domain.  We  suppose  that  both  domains  are  smoothly  bounded. 
Let  Q and  R be  two  points  in  the  intersection  and  let  Q^[u,v]  denote  the 
bilinear  Girichlet  integral  for  the  differential  equation  (2.1.1)  extended 
over  the  domain  D„  It  is  then  easy  to  verify  the  identity 

(3.6.1)  N*(Q,R)- N(Q.R)  » QD+D#[N7r(P,Q)-N(P,Q),N*(P,R)-N(P,R)]  , 
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which  is  a consequence  of  Green's  theorem.  Exactly  the  same  identity,  but 
with  opposite  signs,  is  also  true  if  the  Neumann's  functions  are  replaced  by 
the  corresponding  Green's  functions. 

* 

Let  us  suppose  now  that  the  domains  D and  D are  very  near  to  each  other 
and  that  their  boundary  surfaces  lie  in  an  S neighborhood  of  each  other.  In 
this  case,  in  each  closed  subdomain  of  the  intersection  we  have  a uniform 
estimate 


(3.6.2.)  N*(P,Q)-N(P,Q)  » 0(£  ) 

Hence,  the  contribution  of  to  the  above  Dinchlet  integral  is  only  of  the 

2 

order  £ \ The  first  order  terms  of  the  Dirichlet  integral  are  due  to  the 
difference  shells  and  lead  to  the  equation 

(3.6-3)  N*(Q,R)-N(Q,R)  " J^VN(P,Q)*  VN(P,R)  + p(P)NtP,Q)N(P,R) ] 8iyl<rp+  oU)  , 

C 

?/hich  is  Hadamard's  variational  formula  (2.10.15). 

If  we  want  to  evaluate  the  identity  (6.l)  to  a higher  order  of  precision 
in  £ , we  have  to  utilize  the  following  result.  Let  C be  a sufficiently  smooth 
surface  and  a thin  shell  over  C with  variable  width  Si/.  Let  f(P)  be  a 

twice  continuously  differentiable  function  in  the  closure  ofIZj  then  we  have 
the  estimate 


where  and  jO^  are  again  the  principal  radii  of  curvature  at  the  point  P 
of  integration.  Applying  this  general  formula  to  estimate  (6.l),  we  find 


easily 

(3-6.5) 


O N(QjR)  - 2^1!  8N(P,Q),  $N(P,R) 


£■  rvs>(p 


ex'] 

o ^ / 


vn(p.r)  1 


H(P,Q)N(P,R)~ 


(~jr+  VN(p.Q)‘  VN(P,R)  + fet(p,q)m(p5r) 

r i r 2 J F * 
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It  is  worth  while  to  consider  in  more  detail  the  term  -“  “(VN(P,Q)'  VN(P,R)) 

9 tf 

Let  j(usv)  denote  the  vector  varying  on  the  surface  C and  put 

(3.6.6)  VN(P,Q)  » > VN(P,R)  * Slly 

This  representation  is  possible,  since  VN  has  no  normal  component  on  C.  We 
can  write  ^ 

(3.6.7)  5^rrVN(P,Q )•  V»(P,R)1  - 71  Vi 

i?k»l  I k k 

3 2 

. 7^  a «(p.R)  W,.s).  ,,, 

u£i*%t3\  0Ik  1 ’ 


where  V is  the  normal  vector  to  C, 

(3.6.8)  ^ 9 jj  (£uk  V » 


* r x r 1 
u v1 


We  start  now  from  the  identity 

3 


(3.6.9) 


y7  „ 

ih  3xi  1 


and  we  obtain  by  differentiation 

3 

C 


(3.6.10) 


?2N 


ilSl  0xi  3\  3u 


7S. 


T*  J2JL 


•a*', 

. . . An 

0n 


and  an  analogous  equation  with  v replacing  u.  Thus,  using  (6.6),  we  can  write 

(3.6.111  £ #<“  - - £ eg*  ( J J?) 

k i * 1 i 


. , , 9x,  <9x. 
iak“l  i k 


- - ($  £ + $ r )*()fv  + S 2^  ) 

u r-v  u 4-u  “V 


It  is  well  known  that  the  relations 

(3.6.12)  £ -72  * ~ <£  , r • ~v  = r *2/*  = - 77(  „ £ • -jy  •-  -7? 

u u ’ ~V  U u v v v *- 

hold,  when  ^ 7^,71  denote  again  the  coefficients  of  the  second  fundamental 

form.  Thus,  we  can  put 
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(3.6.13)  ^[VN(P,Q)«VN(P,Rn  = 2 l<Xtftf+Yk(&8  ] 

= 2 7n(p,q)'  vn(p,r)  LioS-  * + $JDl 

v tHJ  v ’ ' 0(tf  E - («  6 ~ )F  + /3  8 G 

Thus,  the  normal  derivative  of  the  scalar  product  between  two  Neumann's  functions 

equals  the  scalar  product  itself,  multiplied  bv  the  ratio  of  the  second  and  the 

first  bilinear  fundamental  forms.  It  is  remarkable  that  his  expression  involves 

only  first  order  derivatives  of  the  Neumann's  functions. 

We  want  to  formulate  our  result  for  the  particular  case  of  Laplace' 3 

equation.  We  find 

(3.6.14)  52N(Q,R)  = aOpCSN^Q^SNCP,!*)] 

+ JJVN(P,Q).  VN(PJR)[2^|VN(P,Q),VN(P,R)|  - (~"  + *jj-)  d<Tp  . 


Here,  “^(x^>Vp)  denotes  the  ratio  of  the  second  and  the  first  fundamental  forms 
for  any  two  given  tangential  vectors  and  y2>  This  result  can  only  be 
applied  in  the  case  that  the  domain  D considered  is  the  exterior  of  a closed 
surface  C,  since  only  in  this  case  can  the  existence  of  a Neumann's  function 
in  the  proper  sense  oe  assumed. 

We  proceeded  in  this  section  in  a rather  formal  way,  without  specifying 
the  exact  assumptions  under  which  the  above  estimates  can  be  justified.  It  is s 
however , sufficient  to  point  out  such  a formal  method  which  works  at  least  in 
the  case  of  analytic  surfaces  C;  for,  we  know  by  our  general  theory  of  variation 
fields  that  a formula  of  the  type  (6.5)  for  the  second  variation  of  the  Neumann's 
function,  must  exist  and  can  be  obtained  by  integration  by  parts  of  the  interior 
variational  formula.  The  method  of  the  present  section  serves  only  to  give  a 
short  cut  to  the’ final  result  without  too  much  labor. 
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Next,  let  ^(P)  denote  the  velocity  potential  for  an  incompressible  and 
irrotational  fluid  flow  around  a surface  C which  has  at  infinity  unit  velocity 
in  the  direction  of  the  x^-axis,  This  function  has  been  discussed  in  Section 
12  of  Chapter  II,  and  its  first  order  variation  was  determined  there.  Let 
C^n  denote  the  corresponding  coefficient  defined  by  (2.12.17),  which  is 

closely  related  to  the  virtual  mass  of  C with  respect  to  the  direction  x^. 

* • 
Let  D be  the  exterior  of  C and  let  C be  another  closed  surface  with  exterior  L 

* * 

such  that  D^*D»D  is  not  empty.  If  (X  ^ and  denote  the  corresponding 
quantities  for  the  new  domain,  we  find 


(3.6.15)  ^11~  ^11  * 477  ^D+D*^ *Pl  ~ ^Pl^  * 

where  Q is  the  Dirichlet  integral  for  Lcplaca's  equation.  This  formula  is 
analogous  to  (6„l)  and  can  be  derived  also  by  application  of  Green's  identity. 
It  is  also  to  be  understood  here  that  and  Cp^  are  defined  to  be  zero  out- 
side of  their  domains  of  definition. 

Evaluating  (6.15)  the  same  way  as  we  evaluated  (6.l),  we  find 


(3.6.16)  S<*u  =»  JJ  d CT 

G 

which  is  already  given  In  (2.12.21),  and 


(3.6.17) 


11 


2tF  Vs  Ti1* 


1 

4TT 


Finally,  we  want  to  give  the  expression  for  the  second  variation  of  „ 
in  the  case  of  an  axially  symmetric  surface  C,  under  the  assumption  that  the 
variation  Si/  is  also  performed  in  an  axially  symmetric  way.  Let  x.y  be  the 
coordinates  in  the  meridian  plane  and  let  y*  0 be  the  axis  of  symmetry.  We 
denote  the  meridian  curve  of  G in  this  plane  by  c.  If  K is  the  curvature 
of  e at  a point  (x,y),  we  have 
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(3.6.18) 


-i.  + _1_ 

A A 


( * + 


1 


) 


for  the  corresponding  points  P on  C,  since  the  domain  D considered  lies  out- 
side of  C.  Using  the  fact  that  because  of  the  axial  symmetry  the  velocity 
vector  V lies  always  in  the  meridian  plane,  we  calculate  easily  by  the 
method  of  this  section  that 


(3.6.19)  (VCp-^2  - "2  KCV^)2 


Hence,  we  have 

(3.6.20) 


s-2  _ 

O (X. 


11  2Tf 


(V<Pi)2(^=i=p--Ky)Sy2as  . 


Consider  finally  the  case  where  C is  a rotationally  symmetric  vortex 
sheet.  In  this  case,  (VCp^)  will  have  equal  values  on  the  upper  and  the 
lower  side  of  the  sheet.  The  curvatures  of  C will  have  to  be  taken  with 
opposite  signs  on  the  upper  and  lower  side*  however,  so  that  in  this  case 
formula  (6.20)  reduces  to 


(3.6.21)  $2fXn  - ^ QdCS  Cpil  >° 

Since,  by  its  very  definition,  the  vortex  sheet  satisfies  the  condition 
8 O(^“0,  we  recognize  that  it  represents  a local  minimum  for  the  coeffi- 
cient 

We  shall  give  another  formula  (4-3-1)  for  the  second  variation  o.f  the 
virtual  ma3s  of  a body  of  rotation.  The  stream  function  4^  will  enter 
instead  of  the  velocity  potential  used  in  (6.20).  Both  formulas  are, 
of  course,  equivalent,  but  (4.3  l)  is  better  adapted  to  the  uniqueness  question 
treated  in  Section  3 of  Chapter  IV. 
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CHAPTER  IV 

AXIALLY  SYMMETRIC  VORTEX  SHEETS 


1. 


A3  an  application  of  the  methods  of  Chapter  II  and  as  a first  step  in 
the  solution  of  the  problem  of  3 -dimensional  vortex  sheets  discussed  in 
Chapter  I,  we  formulate  and  solve  a minimum  problem  which  gives  a construction 
for  vortex  sheets  in  axially  symmetric  steady  irrotational  flow  of  an  incom- 
pressible fluid. 

Let  B denote  a compact  connected  set  in  the  half- plane  y ^0  which 
intersects  the  x -axis,  let  b denote  a point  in  the  infinite  exterior  compo- 
nent of  B relative  to  y > 0 and  let  W be  a continuum  in  the  half- plane  y ^0 
which,  together  with  the  x-axis,  joins  b to  Bin  the  sense  that  the  x -axis  and 
B+W+b  form  a connected  set.  We  denote  by  D the  infinite  region  complementary 
to  B+W+b  in  the  naif- plane  y$s-0  and  we  consider  an  axially  symmetric  flow 
parallel  to  the  x -axis  in  the  3 -dimensional  region  obtained  by  revolving  B 
about  the  x -axis.  This  flow  is  governed  by  a velocity  potential  cp  and  a 
stream  function  ''P  which  are  functions  of  x and  y only  and  satisfy  the 
generalized  Cauchy-  Riemann  equations 


u.i.i)  9X  - ■ %--V* 

in  D and  the  boundary  conditions 

(4-1.2) 

on  the  curves  C bounding  D.  We  normalize  at  infinity  to  have  the  behavior 

(4.1.3) 


dj  2®  „ 

+ ■ ~dv  ' 0 


2 2 

Y „ X-  + 

2 3 
r 


2 2 2 
, r = x 4 y 


and  we  call  C*  the  virtual  mass  of  the  flow. 
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We  consider  B and  b to  be  fixed,  and  we  attempt  to  choose  the  variable 
continuum  W so  that 
(4-1.4)  Of  a minimum 

We  shall  prove  that  a unique  extremal  continuum  W for  (1.4)  exists,  is  an 

analytic  curve  connecting  b with  B or  with  the  x-axis,  and  generates  a surface 

of  revolution  about  the  x-axis  characterized  as  a vortex  sheet  for  the  flow  4* . 

As  a preliminary,  we  point  out  that  only  a bounded  class  of  continua  W 

need  be  considered  for  the  extremal  problem  (4. 1.4).  For  if  W rises  to  a 

height  h above  the  x-axis,  we  can  show  by  symmetrization  of  B+W+b  in  the 

y-axis,  which  decreases  of,  that  the  virtual  mass  pf  of  B+W+b  is  larger  than 

the  virtual  mass  Pf(h)  of  the  flow  pa3t  a disc  of  radius  h perpendicular  to 

the  x-axis  [7,  91.  Since  pf(h)— >00  as  h->oo,  there  is  no  loss  of  generality 

in  requiring  that  B+W+b  lie  in  a suitable  high  strip  0^y*s£:H.  Suppose  now 

that  I is  a curve  lying  within  a very  large  strip  fxl^I,  but  in  no  smaller 

vertical  strip.  Let  4^  be  the  stream  function  of  the  flow  of  the  type  (1.3) 

in  the  exterior  of  the  rectangle  |x|^I,  O^y^H  in  the  upper  half-plane; 

** 

and  let  ^ be  a solution  of  (1 .1),  defined  in  one  of  the  two  rectangles 

obtained  from  |x|^;I,  0^-y^H  by  drawing  a vertical  line  through  b,  which 

has  a suitable  large  positive  value  on  this  line  and  which  vanishes  on  the 

remainder  of  the  boundary  of  the  rectangle.  If  W touches  the  lines  | x|  ■ I 

at  a point  z in  the  half- plane  y>  0 near  z let  'f,+  denote  the  branch  of 
o o 

which  is  defined  outside  the  rectangle  | x|  I,  0^  y^.H,  and  let  V~  be  the 
branch  of  'f*  defined  on  the  other  side  of  W.  Then  ¥ 4^  and  H'”*  4J 

near  zq,  by  the  maximum  principle  for  (l.l).  Hence  at  zq,  for  large  I, 

stsL  > a£.  >^£1  > 

d'U  2>-V  *3  V By 


(4-1.5) 
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since  as  I -^>00  we  have  4*  —^O  and  'r  0 near  z • But  Hadamarcrs  formula 

o 


(4 -1.6) 


St**  2 


(Ji4i)2  -S-iZ-da 
' $2/  7 


can  be  applied  to  show  that  because  of  this  inequality,  a shift  of  the  curve 

+ 

W inwards,  defined  by  bringing  the  vertical  lines  x*-I  closer  together  and 
by  replacing  arcs  of  W beyond  these  lines  by  vertical  segments,  will  diminish  Q\ 
Indeed,  pf  can  be  diminished  thus  until  I is  so  small  that 

(4.1.7) 


aiL  „ 2± 

'dV  dv 


** 


and  hence  we  obtain  bounds  on  the  height  and  width  of  competing  continua  W. 

We  shall  find  it  more  convenient  to  replace  in  D the  stream  function 
by  the  positive  solution 

(4.1.8)  u - ^/y1^2 


of  the  self-adjoint  partial  differential  equation 

(4.1.9)  V2u-(3/4y2)u 

From  a minimal  sequence  of  domains  for  (1.4)  whose  virtual  mass 

coefficients  pfn  approach  their  greatest  lower  bound,  we  find  that  we  can 

select  a subsequence,  again  denoted  by  D^,  which  converges  in  the  sense  of 

Caratheodory  [3]  to  a limit  domain  D.  The  corresponding  functions  u^  form 

an  equicontinuous  family  in  each  closed  subdomain  of  D,  whence  by  suitable 

extraction  of  a further  subsequence  we  can  achieve  that  tend  to  a limit 

solution  u of  (l„9).  Each  u^  is  positive  in  and  hence  is  subharmonic  there, 

by  (1.9).  Thus  if  we.  let  be  a harmonic  function,  defined  in  the  region 

obtained  by  excluding  the  exterior  of  a large  circle  R from  D , which  vanishes 

on  the  boundary  of  D and  is  equal  on  R to  the  maximum  of  u on  R,  then  u ^.U 

n ^ n 7 n n 

* 

in  D . As  the  domains  D converge  to  D,  it  is  known  that  the  harmonic  functions 
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Uq  approach  a harmonic  limit  U which  vanishes  on  the  boundary  of  D„  Hence 
the  limit  solution  u of  (1.9)  in  D must  vanish  on  the  boundary  C of  D,  since 
O^u-^U  in  D*. 

*i  /o 

We  conclude  that  ^ • uy  ' is  the  stream  function  of  a parallel  flow  (1.3) 
in  D,  and  therefore  the  virtual  mass  Q\  for  D lias  the  desired  extremal 
property  (l.4). 

In  the  neighborhood  of  a point  zq  of  the  extremal  continuum  W bounding  D, 
chosen  in  the  half- plane  y > 0 and  not  in  B+b,  we  make  an  admissable  variation 
of  the  form 


(4»1.10)  z » z + , 

z - t 

where  co  is  an  infinitely  differentiable  function  which  is  identically  1 
in  a neighborhood  Cl  of  zq  and  vanishes  identically  outside  a slightly  larger 
neighborhood.  By  Chapter  II,  formula  (2.9.4),  the  variation  of  virtual  mas3  o< 
under  such  a shift  is  found  to  be 


(4.1.11)  $<X  - Re'j  2£  f u2  [“t2 Idxdy-  8rrL  S(t)[-|f]2+  £ Q (t)  f 

L iJa  9z  4y2(z-t)  0 J 

+ o(lel2)  , 

where  QQ(t)  is  analytic  for  t in  O.  and  where  £(t)  is  1 for  t in  D and  0 

for  t exterior  to  D.  Two  applications  of  Green *s  theorem  yield 

r f fP  - I 3u  u.  rt2  3uu  0 

(4.1.12)  Re-}&  JJ  + ■2iL7  + — ^Idxdy- 8TT&£(t)u2+  8 Q(t)| 

DQ  z~  y I6y4  iy' 


+ o(l£l2)  , 

where  Q(t)  is  another  analytic  function  of  t in  Q. 

From  the  extremal  property  (1.4)  of  <X  , we  conclude  that  0 for  all 

sufficiently  small  complex  £.  Hence,  by  the  usual  argument 

~ P T “ 7 3u  u._  n 2 3uu 

(4.1.13)  8i7^(t)ut  » C — g ^ + ^4  + — ^]dxdy+Q(t) 

iSa  y 16f  ^ 
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inO.  We  conclude  from  (1.13 ) that  u2  is  continuous  throughout  0. , if  we 

2 

define  it  to  be  zero  in  the  exterior  of  D,  The  continuity  of  u in  the 

% 

neighborhood  of  each  point  z^  of  W is  the  variational  condition  imposed  upon  W 
by  (1.4)  which  will  lead  us  ultimately  to  deduce  that  revolving  W about  the 
x -axis  generates  a vortex  sheet. 


2. 


2 

Since  u^  is  continuous  inO,  we  can  estimate  the  integrand  in  (1.13) 

2 

once  again  and  prove  that  u^  actually  satisfies  a Lipschitz  condition  there. 
Hence  in  the  region  lu^f  0 we  can  solve  the  differential  equation 

(4.2.1)  u^dt  + " 0 


for  the  level  curves  of  u.  In  particular,  any  arcs  of  the  boundary  u3  0 of  D 
in  the  region  | | must  be  differentiable  curves.  If  we  pick  Cl  to  lie  in 

p 

this  region,  we  can  apply  (1.13)  repeatedly  to  show  that  u*  has  continuous 
derivatives  of  several  orders . Thus  for  a fixed  determination  of  the  square 
roots  involved,  we  show  using  the  proof  of  the  reflection  principle  that 


(4.2.2) 


u 


(u2)l/2  dz+  (u?)1//2  dz^l 

Z Z I 


is  an  analytic  solution  of  (l.9)  in  the  region  I u^  | >0,  and  the  boundary 
arcs  u«0  of  D there  must  be  analytic.  Formula  (2.2)  gives  an  analytic 
continuation  of  u across  these  arcs  according  to  which  u transforms  into  its 
negative  as  we  cross  the  boundary,  and  is  therefore  two- valued. 

In  order  to  see  that  there  actually  exist  arcs  of  W in  the  region  | u , ] ;>  0, 
it  is  sufficient  to  notice  that  every  point  zq  of  W which  can  be  touched  by  a 
small  circle  R lying  in  E is  of  this  type.  For  if  g is  the  Green's  function 
of  (1.9)  in  R with  its  infinity  at  the  center  of  R,  then  for  a small  enough 
£•>0,  we  find  that  in  a neighborhood  of  zq  in  R 

(4.2.3)  u 5;  eg 


Since  u and  g both  vanish  at  zq,  we  can  take  normal  derivatives  there  to 
obtain 


(4.2.4) 


-2ul 

9-v  'dv 


whence  L|gtl'>  0 at  zq. 

It  remains  to  discuss  those  portions  of  W on  which  |u  | *0.  These 
are  identical  with  the  subsets  of  W on  which  |<^’x+i  | “ and  r,e  shall  use 

properties  of  the  quasi- conformal  mapping 


(4.2.5)  » - % + 1 

to  discuss  their  nature.  Let  us  denote  by  P any  region  complementary  to  B+b 
in  which 

(4.2.6)  0<\a>+iV\<£ 

* X X 

for  some  small  E^O.  The  boundary  | (|>x  + i [ * €.  of  this  region  consists 
of  a finite  number  of  analytic  arcs  which  transform  by  (2.5)  into  the 
circumference  |w|*  E , covered  finitely  many  times,  say,  m times.  Since  the 
mapping  w = CD  + i y is  univalent  in  the  small,  we  conclude  that  the  image 

I X X 

of  r in  the  circle  0-c  | w ( <-  6 is  a Riemann  surface  of  m sheets . Hence  the 
boundary  of  P contains  at  most  m continua  on  which  t<px+I(4'x|  * 0,  ami  we 
shall  prove  that  these  continua  reduce  to  isolated  points . 

We  remark  that  the  qua si- conformal  mapping  (2.5)  has  a dilation  quotient 
in  P which  is  smaller  than 

(4.2.7)  K » max  (max  y,  max  ^ ) 

for  all  z»x  + iy  in  P.  If  now  one  of  the  continua  i^x+  i^xi*  0 bounding  P 
were  non- degenerate,  there  would  exist  a positive  lower  bound  S for  the 
^irichlet  integral 


of  any  function  U in  P which  vanishes  on  the  boundary  portions  = ^ 

and.  is  identically  1 on  the  boundary  curves  |(^>x  + ^ e define  such 

a function  U by  the  formula 


r 


0 


(4.28) 


i 


\w  I ^ r 


V. 


log  ^ / log  ^ > 


w I ^ r 


where  We  can  estimate  its  Dirichlet  integral  in  terms  of  the 

quasi-conformal  mapping  w«  <~PX+  iTx  obtoi11  (with  w * u+iv) 

(4-2.9) 


2TT  raK 


J st?  P r ‘f 


, jm-aS- 

(log  *^) 

The  estimate  on  the  right  approaches  zero  as  r — *0,  and  this  gives  a contra- 
diction of  the  hypothesis  that  the  continua  |d|>x'i'  iVx|  * 0 on  the  boundary  of  P 
do  not  degenerate  to  points . 

In  order  to  see  quite  clearly  that  the  continua  | Cb  + iM'*  1 = 0 of  the  set 
71  all  reduce  to  isolated  points,  we  note  that  the  analytic  arc3  of  7!  in  the 
region  |cp  + | ">  0 can  be  given  a certain  order  as  we  proceed  from  b to  B 

along  7t . This  follows  from  the  fact  that  deleting  a small  arc  from  W must 
disconnect  W,  since  it  diminishes  the  virtual  mass  Of.  The  two  disconnected 
pa/ts  of  W axe  composed  of  those  arcs  which  cone  before  and  those  which  come 
after  the  deleted  one,  respectively.  Thus  any  continuum  of  W on  which 

r'41  I * 0 must  be  preceded  and  must  be  followed  by  arcs  of  71  in  the  region 
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|0D  +i^  | and  therefore  must  be  a boundary  component  of  some  region  F"* 

of  the  type  (2.6).  Hence  any  such  continuum  reduces  to  a point,  and  we  can 
conclude  that  these  points  are  removable  singularities  of  the  solution  u of 
(1.9)  defined  by  (2.2).  For  topologically  more  complicated  problems  of  the 
type  (l.4)  than  we  have  posed  here,  the  function  u can  have  a certain  number 
of  branch- points  where  u^  vanishes,  but  in  the  present  case,  because  of  the 
monotonic  dependence  of  (A  on  D,  this  possibility  is  excluded.  The  extremal 
continuum  W must  reduce  to  a simple  arc  without  forks,  since  otherwise  deletion 
of  a branch  of  W,  decreasing  would  be  feasable.  The  arc  W is  an  analytic 
curve,  since  it  consists  merely  of  a level  curve  u*0  of  a regular  solution 
of  (1.9). 

It  follows  from  the  above  arguments  that  the  surface  of  revolution  about 
the  x -axis  generated  by  W is  analytic  and  that  the  speed  IVCp|  0f  the  f low  cp 
past  this  surface  is  continuous  through  it,  so  that  the  surface  represents  a 
vortex  3heet  in  the  flow.  One  checks  that  on  either  side  of  the  vortex  sheet, 
in  the  present  example,  the  velocities  have  opposite  directions.  More  compli- 
cated vortex  sheets  with  several  branches  and  forks  can  be  constructed  by  the 
same  method  with  relatively  little  additonal  difficulty. 


3 . SniqU§R£S£  • 

We  sketch  In  this  section  a proof  of  the  uniqueness  of  the  extremal 

curve  W solving  (1.4)  which  is  based  on  the  ideas  of  Chapter  III.  A.  possibly 

shorter  approach  to  the  whole  problem  might  be  achieved  through  the  Dirichlet 

principle,  but  we  prefer  here  to  emphasize  the  use  of  variational  methods. 

* 

Let  W be  any  curve  joining  b to  the  set  B or  the  x -axis.  The  curves  W 

* 

and  W , together  with  B and  the  x "axis,  bound  certain  subregions  D^,  of  D. 

In  each  region  T>V)  we  let  U be  the  harmonic  function  which  vanishes  on  the 
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arcs  of  W bounding  D^,,  which  has  the  value  1 on  the  arcs  of  W bounding 

and  which  ha3  a vanishing  normal  derivative  at  any  points  of  B or  the  x-axis 

bounding  Dv„  We  denote  by  the  aggregate  of  all  level  curves  U ■ t in  the 

various  regions  Oy  and  we  denote  by  <X  the  virtual  mass  of  the  flow  (1.3) 

* 

past  the  object  B+W^+b,  In  particular.  W *W  and  W^»W  , while 
and  a ^ Is  the  virtual  mass  coefficient  corresponding  to  W . 


From  the  formula 
(4.3.1)  S2(X 


/ 


+ jj?  £ S u)^  + Sui^ldTr 


for  the  second  variation  of  $ under  a normal  shift  of  the  boundary, 
we  derive  the  expression 


(4-3.2) 


d2CX\ 


dt* 


du.  o o du  0 

t(V  J*)  + 4^  (-^t)2ld'C 


4y 


2 'dt 


for  the  second  derivative  of  Of^.  By  the  continuity  of  (Vu)  across  Wq  * flf 

and  by  (1.6),  we  have 
dor 


(4.3.3) 

I 

dt 

= 0 

5 

t»0 

and  by  (3.2) 

d CX\ 

(4.3.4) 

, X 

2 

1 > 0 , 

O^t^l 

dt 

It  follows  that  (X, 

t 

>CX  for 

t >0,  and,  in  particular, 

(4.3.5) 

cx*  2 

> Of 

J 

* 

provided  W 

does  not 

, coincide 

with  W.  Formula  (4.3.5) 

extremal  curve  W for  (1.4)  is  unique.  A closer  examination  of  the  proof  even 
shows  that  an  arbitrary  axially  symmetric  vortex  sheet  joining  b to  B has  a 
meridian  curve  satisfying  (l»4),  and  hence  a vortex  sheet  such  as  we  consider 
here  must  be  unique. 
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CHAPTER  V 

VARIATION  OF  EIGENVALUES 


1. 


We  consider  now  the  eigenvalue  problem  of  determining  a non- trivial 
solution  of  the  partial  differential  equation 

(5.1.1)  ^u+  Au  - 0 


irhich  vanishes  on  the  boundary  C of  the  plane  domain  D.  We  consider  the 
eigenvalues  A as  functionals  of  the  varying  domain  D. 

Hilbert  [11]  proved  the  continuous  dependence  of  the  eigenvalues  A on 
the  coefficients  of  a differential  equation  if  the  basic  domain  D is  kept  fixed. 
We  can  utilize  his  method  and  his  result  by  the  artifice  used  in  Chapter  II. 

We  transform  the  domain  D into  a domain  D by  msans  of  a deformation  (2.2.2) 
and  ask  for  the  eigenfunctions  u (x^,Xg)  an<*  eigenvalues  A of  the  same 
equation  (l.l)  with  respect  to  the  new  domain  D . Then,  we  introduce  the 
functions 

(5.1.2)  UCx^^jS)  « u*(x*,x*) 

defined  in  the  fixed  domain  D and  treat  the  eigenvalue  problem  for  this 
fixed  domain  arising  from  the  transformation.  We  may  describe  U(x  as 

an  eigenfunction  in  D with  respect  to  the  differential  equation  (see  (2.2.9) 
and  take  p*»yD  « 0) 

(5.1.3)  L.[U]  + A*0U  - 0 

£ 

and  the  boundary  condition  U-0  on  C.  In  this  way,  the  domain  dependence  of 
the  eigenvalue  A is  translated  into  its  dependence  on  the  coefficients 
of  the  new  differential  operator  L^. 
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We  may  now  proceed  as  in  Chapter  II  and  derive  the  variational  formulas 
for  the  eigenvalues  A by  a proper  use  of  Green's  Identity.  We  may  assert 
in  view  of  (2.2.8)  that  all  coefficients  of  the  differential  expression 
depend  analytically  upon  the  parameter  £ . It  is  then  easy  to  see 
that  all  non-degenerate  eigenvalues  A and  their  corresponding  eigenfunctions 
U depend  analytically  on  £,  provided  U is  normalized. 

In  order  to  prove  this  statement  we  observe  at  first  that  the  Green's  .function 
g(p>Qi  £ ) belonging  to  the  differential  operator  1.^  depends  analytically  on  £ . 

In  fact,  the  parametrix  function  s(P,Qj&)  of  this  differential  operator, 
defined  for  two  independent  variables  in  analogy  with  (2.3.5),  may  be  chosen 
to  depend  analytically  on  8.  for  P/Q.  The  Green's  function  g(P,Q; £ ) may  then 
be  obtained  in  form  of  a Neumann's  aeries  which  solves  the  integral  equation 
(2.4.20)  with  the  analytic  kernel  (2.4.21)  in  6 . Since  the  Neumann's  series 
converges  uniformly,  we  see  that  g(P,Qjg)  depends  analytically  upon  this 
parameter  also,  at  least  in  a neighborhood  of  the  value  (L  • 0. 

The  eigenfunctions  U(P;£)  of  the  differential  equation  (1.3)  may  also 
be  considered  as  eigenfunctions  of  the  integral  equation 

(5.1.4)  U(p5e)-  A*jj  g(p,Qje)u(Qje)e(Q)dd:Q  - 0 

D 

# ■ 

■with  tbs  same  eigenvalues  A • The  eigenvalues  A appear  as  the  roots  of 
an  entire  function  in  A which  is  also  analytic  for  sufficiently 

small  £..  If  A is  non- degenerate,  we  have  0 near  the 

value  £ » 0 and  A * i ( £ ) is  an  analytic  function  of  £.  Since  the 
eigenfunctions  U(Pj  £, ) can  be  represented  as  Fredholm  minors  which  depend 
analytically  on  £ , we  have  thus  also  proved  the  analytic  dependence  of 
U(P;£)  on  £ , if  the  corresponding  eigenvalue  is  non-degenerate. 
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We  will  apply  this  result  in  particular  to  the  case  of  the  lowest  eigen- 
value of  the  partial  differential  equations  (1.3)  and  (l.l)  and  its  corresponding 
eigenfunction.  It  is  well  known  that  the  lowest  eigenvalue  of  an  equation  (l.l) 
is  non-degenerate  ana  consequently  we  may  develop  U(P$  £ ) and  A ( S ) into 

V*,  „ ^ 

power  series  in  E which  converge  for  small  values  of  £ . Henceforth,  A \>t,j 
and  U(P;6)  shall  always  denote  the  lowest  eigenvalue  and  its  corresponding 
eigenfunction.  It  Is  obvious  how  to  generalize  the  following  considerations 
to  the  case  of  a non- degenerate  eigenvalue  of  arbitrary  order. 

Since  we  can  put 


(5.1.5)  A*(e)  - + l E2  A2+ ... 


(5.1.6)  U(PjE)  - uo(P)  + Eu1(P)+^  E^u^P)  + ... 

where  A^  and  uq(P)  are  the  eigenvalue  and  eigenfunction  of  the  original 
domain  D,  we  may  calculate  all  other  terms  A^  and  u^/P)  of  the  above  series 
development  by  inserting  these  series  into  (1.3)  and  comparing  the  coefficients 
of  equal  powers  of  £ . This  procedure  is  called  In  physics  the  "perturbation 
method"  and  is  widely  used  in  applications . It  is  easy  to  handle  and  satis- 
factory in  most  problems  in  applied  mathematics.  Obviously,  it  works  only  in 
the  case  of  analytic  coefficients  in  the  differential  equation.  We  can  now 
apply  this  procedure  in  order  to  study  the  dependence  of  the  eigenvalues  upon 
their  domain  of  definition. 

We  introduce  for  this  purpose  the  tensor 

3u  d u 


(5.1.7) 


Vp) 


2 


which  is  symmetric  and  satisfies  the  equations 

(r.i.8) 

k-1  0 axi 


U 
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We  obtain  the  following  formula  similar  to  (2.6,6) s 


It  is  easily  seen  that  by  virtue  of  the  equations  (1.8)  the  integrand  is  a 
divergence  term  and  that  the  expression  (1.9)  may  be  reduced  to  an  integral 
extended  over  the  boundary  C of  D.  In  this  way,  we  obtain  in  analogy  with 
(2.6.15) 

(5.1.10)  SA  * \ ds 

C 


This  formula  is  due  to  Rayleigh  and  is  derived  here  from  the  method  of 
interior  variations.  If  we  use  the  complex  notation  (2.8.6)  introduced 
in  Section  8 of  Chapter  II,  we  find 


(5.1.11) 


d A 
d£ 


- Re 


£-0 


+ 2 ^ ( A Fjujhd'C 

^7  Z O O 


Let  us  remark  finally  that  the  formula  (l .10)  exhibits  clearly  the 
monotonic  dependence  of  the  eigenvalue  A upon  the  domain  D.  Vfe  see  that 
for  8v  >0,  that  is,  for  a shrinking  domain  D,  the  eigenvalue  A increases. 
Analogous  formulas  can  be  given  for  all  higher  non- degenerate  eigenvalues  of 
the  differential  equation  (3„l)  and  similar  conclusions  can  be  drawn. 


2 • The -Hilbert-Green  function . 

In  this  section  we  will  calculate  the  first  coefficient  u^(P)  in  the 
development  (1.6)  of  the  eigenfunction  U(P;£).  In  this  way,  we  will  also 
prepare  the  determination  of  the  term  A 2 which  will  give  the  formula  for 
the  second  variation  of  the  eigenvalue. 

In  order  to  calculate  u^,  we  define  the  Hilbert-Green  function  P(P,Q) 
for  the  differential  equation  (l.l)  and  the  domain  D as  follows s 
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a)  The  function  r(P»Q)  Is  twice  continuously  differentiable  In  B, 
except  at  P*»Q.  It  satisfies  the  inhomogeneous  partial  differential  equation 


(5.2.1) 


b) 


in  D . 


L[n(p,Q)]+  A P(p,q)  - vfr(p,Q)+  A p(p,q) 

O 0*0 

*“  u (p)u  (q) 
o o 

The  function  + log  PQ  is  continuously  differentiable 


c)  P(P,Q)  vanishes  for  P£C,  Q€D„ 

The  Hilbert-Green  function  has  to  be  applied  to  differential  operators 
for  which  a proper  Green's  function  does  not  exist,  as  is  the  case  for  (l.l). 
The  existence  of  this  function  can  easily  be  3hown  from  the  basic  theorems 
on  integral  equations.  In  fact,  let  G(P,Q)  be  the  Green's  function  of  D 
with  respect  to  the  differential  equation  L^Cul^O.  We  have  obviously 


(5.2.2) 


u (P) 
o 


G(P,Q)u)(Q)d'CQ 


9 


and  r(P,Q)-G(P,Q)  « H(P,Q)  is  a continuously  differentiable  function  in  D 
which  vanishes  for  P6C.  H(P,Q)  satisfies  the  differential  equation 


(5.2.3)  75(P,Q)+  A H(P,Q)  - u (P)u  (Q)-  A G(P,Q) 

r o o o o 

If  we  can  show  conversely  that  (2.3)  has  a continuous  solution  H in  D which 
vanishes  for  Pe  G,  we  can  construct  P(P,Q) * G(P,Q) + H(P,Q) . The  necessary 
and  sufficient  condition  for  the  existence  of  a solution  of  the  inhomogeneous 
differential  equation  (2.3)  which  vanishes  on  the  boundary  C of  B is  the 
orthogonality  of  its  right-hand  side  to  the  eigenfunction  uq(P).  But  this 
condition  is  fulfilled  in  view  of  (2.2)  and  the  existence  of  the  Hilbert- 
Green  function  has  thus  been  established. 
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We  observe  that  the  above  conditions  (a)j  (b),  (c)  do  not  yet  determine 
r(P>Q)  uniquely,  since  the  addition  of  any  multiple  of  uq(P)  to  F(P,Q) 
does  not  affect  these  requirements.  We  determine  P(P,Q)  uniquely  by  the 
additional  condition 

(5.2.4)  n p(p  ,Q)uo(Q)dtTQ  ■ 0 

D 

It  can  be  shown  in  the  usual  way  that  the  Hilbert-Green  function  P(P,Q)<, 
which  is  now  unique , is  symmetric  in  both  Its  argument  points . 

Consider  now  a function  v(P)  which  satisfies  the  inhomogeneous  equation 


(5.2.5) 


L [v]  + A v(P)  - f(P) 
0 o 


'(p)  « 0 for  Pec 


We  have  by  Green's  identity  and  the  vanishing  of  P and  v on  C 


(5.2.6) 


1 


vLo[P]-  PLo[v]JdrQ  - v(P ) 


Using  now  (2.l)  and  (2.5),  we  obtain  the  following  representation  for  v(P): 


(5.2.7)  v(P ) - -fj  r(p 

This  result  is  analogous  to  the  solution  (2.1.6)  of  the  corresponding 
inhomogeneous  equation  (2.1.5)  in  the  case  when  a proper  Green's  function  exists. 
We  observe  that  equation  (2.5)  determines  v(P)  only  up  to  a multiple  of  the 
eigenfunction  uq(P),  which  is  clearly  shown  by  the  solution  formula  (2.7). 

We  remark  also  that  f(P)  cannot  be  prescribed  arbitrarily,  but  must  be 
orthogonal  to  the  eigenfunction  uq(P). 

For  later  applications  we  want  to  prove  that  the  quadratic  functional 


,Q)f(Q)d/r0  + u (P) 

H o 


ff- 


(QMQ)d'tr 


Q 


(5.2.8)  -JJ  JJ  P(P,Q)f(P)f(Q)d'UpdarQ 

D D 


based  on  the  Hilbert-Green  function  is  positive  semi-definite  and  that  it 
vanishes  only  if  f(P)  is  a multiple  of  the  eigenfunction  u^(P).  We  may 
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put  every  function  f(P)  into  the  form 


(5.2.9) 


f(P)  - uo(P)J  J fudr  + fx(P) 
D 


where  (P)  is  orthogonal  to  u^CP).  Because  of  (2.4),  we  have 

(5.2.10)  - /-j 

so  that  our  assertion  will  be  proved  if  we  show  that^^fjfJ-  is  positive- 
definite  for  all  functions  f(P)  which  are  orthogonal  to  uq(P). 

For  all  such  functions  f(P)  the  inhomogeneous  equation  (2.5)  has 
solutions  v(P)  of  the  form  (2.7).  Hence,  we  may  write 

(5.2.11)  * “Jj P(p)v(P)dtTp  - - fJ(LoCv]+  AQv)vd'C 

D L 

Applying  finally  Green's  identity,  we  obtain 

(5.2.12)  /-[f,fj  - JJ[(Vv)2-  A^ld'C 

It  is  well  known  that  we  can  characterize  the  first  eigenvalue  of  the 
differential  equation  (l.l)  a3  the  minimum  value  of  the  ratio 

2. 


II 


dTT 


for  all  continuously  differentiable  functions  v(P)  in  D which  vanish  on 

the  boundary  C.  The  minimum  value  is  accepted  only  if  v(P)  is  a multiple 

of  the  eigenfunction  u (P).  Hence,  we  conclude  from  (2.12)  that 

o 

(5.2.13)  «/£f,f^ 

for  all  functions  f(P)  which  are  orthogonal  to  uq(P).  Equality  can  hold  in 

(2.13)  only  if  v(P)  is  a multiple  of  uo(P),  which  implies  by  (2.5)  that 
f(P)  HS-0 . This  shows  the  po3itive-def inite  character  of^^f,^-  . 
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We  apply  now  the  solution  formula  (2.7)  to  the  equation 

?S  7>S 

(5.2.14)  L [il]+  A u _ * - Anu  - l/[u  ]••  A u (~r— ^ 

oT.  ol  looo  o 0 

obtained  by  comparing  £•  terms  in  the  relation  found  when  we  insert  (1.5) 
and  (l.6)  into  (l.3)>  This  g3ves 

^ (i  ^ 

(5.2.15)  u^P)  » a^Cuo(Q)]+  ^oUo(Q)(a|J  + FJ^)]  rXP,Q)d‘T7Q+  JJlu.o(P) 


where 

(5.2.16) 


D 


0S.  V*. 


We  can  simplify  the  result  (2.15)  if*  we  introduce  the  symmetric 
variational  tensor 


(5.2.17) 


T (P.a) . iv(p)  a Onal . 9u°(F)  arinsa 

vik^’8J  aXl  3^  9^ 

■ sikVuo(p).vr(P,Q) 


We  may  then  bring  (2.15)  into  the  form 


(5.2.18) 


ds 


2 


+ WU  (Q) 

I 0 


We  have  the  differential  equations 
2 


(5-2.19)  XI  - - A.  Sr  tu  (p)  P(p;q)i  + J u (q)  u (p)2 


k*  1 ik 


o 2x^  ' o 


2 o'^'  3xt  o 


Hence,  we  may  bring  (2.18)  into  the  divergence  form 

2 r 2 


(5.2.20)  «,(«)-]]  C/-j  II 

K 1 K 11  1 v J 


1 

2 


p p 2_^ 

u (Q)  \ \ X^  (Su2)dt 
0 J J kTi  ^*k  k ° 
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We  can  reduce  this  domain  integral  to  a line  integral  extended  over  the 
boundary  curve  C of  D.  We  use  the  fact  that  and  u^  vanish  on  C and  take 
into  account  the  singularity  of  the  tensor  V at  the  point  Q.  We  obtain 

3u  (Q)  0 op( p n)  ^ » 

(5.2.21)  UX(Q)  - iLj  8,(0)  - j gv’®  gVp  (S-V)dSp  • 

X — Q 

Let  us  return  to  the  eigenfunction  u (Qj'o)  of  the  original  differential 

* 

equation  (l.l)  with  respect  to  the  variable  domain  D . Since 

(5.2.22)  U(x1,x2;  8.)  3 u*(x1+£S1,x2+£S2;  6)  , 


we  derive  from  (2.2l) 
(5.2.23)  ^ ic*" 


or,  in  equivalent  variational  form 


0u  (P) 

Q 


(S-^)ds, 


(5.2.24) 


8 u(Q) 


Svdi 


The  function 


(5.2.25) 


h(Q) 


3u  (P) 


afer  (3'y)dsp 


will  play  a role  in  the  second  variation  of  the  eigenvalue  A . We  observe 
that  because  of  (l.lO)  and  (2.l) 


(5.2.26)  Loth(Q)]+  Aoh(Q)  - - A^Q) 

and  we  verify  easily  that  h(Q)  has  on  C the  boundary  values 

0u  ^u(Q) 

(5.2.27)  h(Q)=-— (S*y)  - - 2_^  S (Q) 

^ v i-1  1 ^Ji 
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3.  Sflgon^  jtaidajblQn.-Q£  -Slgamlag . 

We  are  now  able  to  calculate  the  coefficient  ^ the  series  development 
(1.5)  for  the  eigenvalue  A ( £ ) under  a deformation  (2.2.2)  corresponding  to 
a normal  shift  of  C,  It  is  possible  to  achieve  by  a laborious  process  of 
rearrangement  of  terms  and  of  integration  by  part3  a rigorous  derivation 
very  similar  to  those  performed  in  Chapter  II.  The  end  result  of  the 
calculations  can  be  obtained  easily,  however,  by  the  following  heuristic 
argument . 

■g. 

We  suppose  that  C is  an  analytic  curve  and  that  u(P)  and  u (P)  are  both 
regular  functions  in  the  original  domain  D and  satisfy  the  differential 
equations 

(5.3.1)  Au+  A u ■ 0 , Aui  + A*u*  - 0 


in  D.  We  apply  Green’s  identity 
(5.3.2) 


jj  |u*[A2u+  A u]-  uC du  + A*u*]J  dr 


(u  ~ u*  -^7)ds+  (A  - K) 


* 

uu  dT 


Since  u(P)  vanishes  on  C and  in  view  of  (3.l),  we  obtain 

(5.3.3)  ( A - A )Jj~  uu  dt  - f u “g-y" 

JD  0 

We  observe  that  because  of  (2.23)  and  (2.25) 

(5.3.4)  u*(P)  - u(P)  + £ h(P)  + o(£  ) 

# * 

Assuming  u (P)  vanishes  on  a deformed  curve  C , we  have 

(5.3.5)  0 - u (P  ) = u (P)  + N(s)  + \ £2k(3)2  + . . . 

where  £N(s)  represents  tue  amount  of  the  normal  shift  from  C to  C at  the 


point  P . Hence 
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(5.3.6)  ( A - A )[l  + £ \ \ hud'll  + o(£ ) ] 

JD 

" f£  ft  ^ l-A^a.ote2) 


c 


\ ( OSL.}*'  p ji/jo  ^ 1 \ -2 — IL  .2  VL  ^ — S-b-  TMg 

" J 13t/]  bMu  2 J 9j/2  ?Z/  b * d“  b 9V  dV  ~s 


+ o ( t2) 


We  observe  that  the  last  right-hand  integral  can  be  written  in  view  of  (2.27) 
and  (2.26)  as 

(5.3.7)  £2  hds  " “ ^[(Vh)2-  A^ld-c  + A1  hudZ 

C D JD 

2 

Finally,  we  put  A » ^+^i  + ^2~  ^2+  and  by  1138  of  ^■Lo^0^  and 

comparing  the  coefficients  of  £ on  both  sides  of  (3.6) 

2 

(5.3.8)  |AN(a)2d3-2^t(Vh)2-  A/HT  . 

C 9 D 

Since  u(P)  satisfies  the  first  equation  (3.l)  and  vanishes  on  C,  we  have 

(5.3.9) 


Jlin 


_ JL  2jl 
2 p dv 


where  p(s)  i3  the  radius  of  curvature  at  the  point  P with  the  parameter 
value  s.  Hence,  we  may  bring  (3.8)  also  into  the  form 

(5.3.10)  A2 -- $<f£02  a^<te-2^c(Vh)2-  K/w 


D 


In  general,  h(F)  will  not  vanish  identically  on  C and  we  cannot  assert  that 
the  second  right-hand  integral  is  non- negative . Thus,  even  in  the  case  of  a 
convex  curve  C (p^-0)  we  cannot  be  sure  that  X2  ^-s  negative. 

We  shall  consider  In  the  next  sections  variations  of  a more  special  form 
for  which  the  convex  dependence  of  X upon  the  variation  parameter  can  be 
3hown.  It  will  appear  that  interior  variational  formulas  are  more  convenient 
for  this  purpose  than  formula  (3.10), 
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4-  C^Jivesiiy. 

The  variational  formulas  for  the  eigenvalue  A and  the  eigenfunction 
u(P)  become  very  simple  if  we  assume  that  the  complex  vector 


F(z,z)  « S^fx^x,,)  + iS^XpX^) 

is  an.  analytic  function  of  the  complex  variable  z.  In  this  case,  we  have  by  (l.ll) 


(5.4.  l) 


lA 

dE 


- -2  Aq  Re  J ffp'(a)uo(js)2d'C 
8«o  Id 


o p # 

The  analogous  formula  for  d A/d  £ at  6*0  under  a variation  z *z+£F(z) 
is  more  difficult  to  derive  from  (3.10),  because  such  a deformation  is  not 
equivalent  to  a normal  shift  when  second  order  terms  are  involved.  But  for 
a conformal  mapping,  Lg.[U]  = and  ©»|l+  £F/12,  so  that  substitution 
of  (Lo5)  and  (1.6)  into  (1.3)  and  comparison  of  terms  in  yields 

(54.2)  V2V  A*,  - -2>>1U1-4A0U1  He  {?'}  - AjU,  - 4 h?0  Re{F']  ' 2 Vo'F' |2 


Since  u,  vanishes  on  C,  the  eigenfunction  uq  must  be  orthogonal  to  the  right- 
hand  side  of  (4.2),  and  this  gives 

(54.3)  A2  - He.{V}  « -2A0^%|r'!2dr  -a^Vi" 

D D D 

UsiDg  the  condition 

(5.4-4)  + A u » - -2/  u Re  |f'"1 

1 0 1 ''1  0 o o t j 

for  £ terms  in  (1.3)  to  cancel  and  noting  that  with 

(5.4-5)  f » ,A,u  +2  A u Rs-Tf-']. 

1 o 00  L J 

foricula  (2.12)  gives 

(54.6)  rflf.f]  - jJuVu/-  A 0uJ]dT 

D 
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we  bring  (5.4.3)  into  the  form 


(5.4.7)  A2  = -2Ao5j^lF'|2dir  -4^^^  Re-jVj  d^T 


Using  (4.l),  we  simplify  (4.7)  to  obtain 


(5.4.8)  X M?  - 2 (#)' 


if 


‘d£ 


£ =0 


-2  A/ - 2 jjlF'l2 


We  observe  that  the  right-hand  side  of  (4.8)  is  always  negative,  except  for 
the  case  F/(z)  = 0;  the  deformation  induced  in  this  exceptional  case  is  a 
translation  of  the  whole  plane,  and  it  is  obvious  that  in  this  case,  indeed, 
A ( £ ) ■ A(0).  Thus,  we  are  led  to  the  differential  inequality 


(5.4.9) 


d V 


[* 


Mil 


1 


[ 


Mr-  2 


A3  ‘ de2  " >de 

whicn  proves  the  convexity  of  the  reciprocal  of  the  eigenvalue  under 
infinitesimal  analytic  deformations. 

In  order  to  generalize  this  result  to  finite  deformations,  we  consider 
the  partial  differential  equation 

(5.4.10)  V2u+  A |h'(z,T)l2u  ■ 0 , u » 0 on  G , 


where  h(z,T)  is  an  analytic  function  of  the  complex  variable  z6D  and  the 

real  parameter  T,  O^T-^1.  We  3et  h'  » and  assume  h/(z,0)»l. 

o z 

Obviously,  (4.10)  is  the  membrane  equation  for  a domain  obtained  from  D 
by  a mapping  h(z,T),  referred  back  to  the  domain  D.  We  observe  that  the 
differential  equation  can  be  treated  independently  of  the  univalency  of 
h(z,T)  and  permits  a definition  of  eigenfunctions  and  eigenvalues  of  the 
membrane  problem  also  for  non-schlicht  domains.  We  normalize  the  eigenfunctions 
u(z;T)  of  (4.10)  by  the  condition 

^^u2|h/  (z,T)  j^dT  - 1 


(5.4.11) 
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and  consider  u^-ufzjT)  and  AT»  A(T)  in  their  dependence  on  T.  We  may 
repeat  all  calculations  of  this  chapter  leading  to  expressions  for  and 

^2  t 

— 2 , but  now  these  derivatives  can  be  determined  for  finite  values  of  T. 

d'T 

We  introduce  a Hilbert-Green  function  PT(z, % ) by  the  conditions 

(5.4.12)  ATlh/(z,T)|2  rT(ztZ)  - lh'(z,T)l2uT(z)uT(^)  , 

(5.4.13)  rT(z,^)  +^iog|z-£i 

is  continuously  differentiable  in  D , 

(5.4.14)  rT(z,^)  - 0 for  zee,  £eD  , 

(5.4.15)  JJ  r^,(z,  ^ )u,j,(z)l hx  (z,T)|  2d^  "0 

D 

The  existence  of  such  a fundamental  function  can  be  shown  in  exactly  the 
same  manner  as  w as  done  in  Section  3.  Let  be  the  quadratic  functional 

defined  in  (2.8),  but  with  instead  of  P as  kernel.  W©  see  immediately 
that  X?  "{f’f}  is  also  semi-definite  and  vanishes  only  for  the  eigenfunction  u,p. 

One  finds  by  easy  computation 

(5.4.16)  - -AjjuT(z)2  ~ {\h/{z,T)\2)A'C 

and 

(5.4.17)  ^2'“  2 ’ X(f)'  ^dT*^  + 2o^T'lfT,fT}  + ^t2  "jjh  (z/1')!  '^(zfdt;  - 0 , 

with 

(5.4.18)  fT  . | h^z.T)!2*  A(T)  -Jj  (I  h/(z,T)  12)1 

Thus,  we  may  write 

(5.4.19)  ^ ^"a}tP  +'A“  jj  ^t2  lh^z»T)'2 
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The  convexity  of  A(?)  is  thus  ensured  if  ne  can  assert  that 

(5o4.20)  |h/(z,T)|2  ^0,  for  zeD  , 

d T 

This  is,  for  example,  the  case  if 


(5.4.21)  h(z,T)  = 2 + 17(2) 


since  we  have  then 
2 

(5.4.22)  ^2  |fc/(z,T)|2  - 2|F/(z)|2 

We  have  shown  that  the  reciprocal  A(T)  ^ of  the  eigenvalue  of  a family 
of  domains  is  a convex  function  of  the  parameter  T if  the  domains  are 
obtained  from  the  original  domain  D by  the  conformal  mapping  (4.2l).  It  is 
not  necessary  that  all  domains  Hj  be  schlicht  in  the  complex  plane. 


5 . Application  to  extremal  problems . 

As  an  application  we  consider  the  family  3“  of  univalent  functions  f(z) 
in  the  unit  circle  with  f(0)»0,  f'(o)»l,  and  aslc  for  a function  £ (z)€d~ 
which  leads  to  a domain  A with  maximum  first  eigenvalue  A of  the  membrane 
problem.  We  introduce  the  function  h(z  ,T)  ■ z * T[f^(z)  - z1-  ar.d  the  corresponding 
eigenvalue  A<j.*  We  find  from  (4.16)  and  (4.19) 

(5.5.1)  £ (|)  .alf'-lfld-C  . 

D ^ J 

(5.5.2)  ^2  (-^-)  - “^2^T  + ‘f'JJuT  Ujj,"1*  dT  ^*0 

Observe  that  the  first  derivative  of  vanishes  for  T«0j  this  follows  from 
the  fact  that  f'Tz)-  1 vanishes  at  the  origin  and  from  the  radial  symmetry 
of  the  first  eigenfunction  for  the  circle.  Hence  by  (5.2),  4?  (l/A  ) will  be 
positive  for  all  T > 0 and  we  find 


y 


(5.5.3) 
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if  A is  the  lowest  eigenvalue  of  the  unit  circle.  But  since  A(l)  is 
the  largest  possible  value  of  A , ?;e  conclude  that  fm(z)=z;  that  is,  among 
all  equivalent  domains  the  unit  circle  has  the  largest  eigenfrequency . This 

II 

result  is  due  to  Polya  - Sze go  [16  ]„ 

We  obtain  the  following  corollary*  if  a domain  D satisfies  the  condition 

(5.5.4)  J|u2F(z)d'T  - 0 

D 

for  all  functions  F(z)  regular  in  D and  vanishing  at  a point  P£D,  where  u 1b 
the  first  eigenfunction  of  D,  then  D is  necessarily  a circle  with  P as  center. 
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